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Executive  Summary 


This  final  report  contains  a  summary  of  the  principle  finding. 
Experimental  and  Theoretical  details  are  presented  in  greater  detail  in  a 
series  of  Appendices.  This  format  was  chosen  since  this  report  may  be  read 
by  individuals  with  specific  interests  thereby  enabling  the  reader  to  find 
the  relevant  detailed  information  in  one  place. 

The  principle  research  objectives  of  the  grant  were: 

1)  Investigate  the  infrared  and  optical  receiving  properties  of  STM 
irradiated  tunnel  junctions  and  to  determine  the  effects  of  geometry, 
temperature,  and  materials  properties  on  the  receiving  properties. 

2)  Investigate  the  proposal  that  the  electron  transit  time  in  irradiated 
vacuum  tunnel  structures  can  be  characterized  in  terms  of  the  natural  time 
scale  provided  by  the  laser  frequency  and  to  determine  the  physical 
mechanisms  limiting  the  response  of  the  tunnel  device. 

Several  of  the  significant  accomplishments  achieved  were: 

1)  Successfully  fabricated  atomically  sharp  silicon  emitter  embedded  in 
a  thin,  30  Angstrom,  tunnel  oxide.  This  represented  the  first  fabrication  of 
the  solid  state  equivalent  of  a  STM  as  a  circuit  element. 

2)  Extended  the  fundamental  understanding  of  the  limitations  of  the 
planar  Fowler-Nordheim  theory  when  applied  to  tunneling  from 
atomically  sharp  emitter  tips  by  calculating  the  fields  induced  at  these  tips 
using  a  more  realistic  prolate-spheroidal  co-ordinate  model.  This  is  a 
marked  deviation  from  the  piurely  planar  Fowler-Nordheim 
approximation.  This  calculation  represented  the  solution  to  a  previously 
unsolved  problem  in  electrostatics.  A  fully  three  dimensional  analytic 
solution  for  a  sharp  tip  and  a  planar  base  was  foimd.  The  analytic  solution 
included  the  effects  of  the  image  and  multiple  image  potentials  in  the 
model  calculation. 
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3)  Calculated  the  effects  of  inelastic  processes  in  time  dependent 
tunneling  in  an  STM  junction  using  a  kinetic  formulation  of  the  transfer 
Hamiltonian  method. 

4)  Designed  and  built  a  digitally  controlled,  network  accessible  STM 
which  affords  remote  accessing  over  the  university  network  and  real  time 
data  acquisition.  All  aspects  of  the  experimental  control  of  the  STM  are 
under  computer  control,  including  coarse  and  fine  approach,  settling  and 
pulse  control  of  laser  and  incoherent  sources. 

5)  Determined  that  the  role  of  incoherent  photoassisted  tunneling  in  the 
visible  region  is  minimal.  Incoherently  stimulated  tunnel  currents  were 
below  0.01  nanoamperes.  This  is  contrasted  with  typical  STM  tunnel 
currents  on  the  order  of  1  nanoampere. 

6)  Preliminary  experiments  of  the  role  of  surface  modification  in  the 
Fowler-Nordheim  or  field  emission  range  of  tunneling  appears  to  be 
associated  with  a  field  effect  and  not  with  thermal  evaporation.  This  has 
important  implications  in  the  fabrication  of  10-20  Angstrom  feature  sized 
devices.  Experiments  are  continuing  on  establishing  the  mechanism  for  the 
fabrication  of  nanosized  dots  on  the  order  of  20  Angstrom  diameters  and  10 
Angstrom  depth. 

7)  Developed  the  imaging  and  control  software  to  monitor  surface 
modification  between  laser  pulses  to  establish  greater  accuracy  in  nano 
distance  measurements  which  affect  determination  of  tunnel  times. 


These  investigation  have  important  payoff  to  the  Air  Force  in  the  areas  of 
surface  modification  and  nano-scale  fabrication  and  the  area  of  high  speed 
tunnel  devices  acting  as  either  switching  elements  or  as  infrared  and  optical 
detectors. 
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Introduction 


Selective  room  temperature  infrared  and  laser  receivers  incorporated  into 
silicon  integrated  circuits  have  important  implications  in  high  speed 
infrared  signal  processing  and  optical  computation.  Although  a  wide  range 
of  ni-V  and  II-AA  compound  semiconductors  exist  with  stoichiometrically 
adjustable  band  gaps  these  materials  must  be  heteroepitaxially  grown  by 
MBE  or  MOCVD,  and  as  such  are  costly  and  difficult  to  produce.  In 
addition  long  wavelength  infrared  detectors  using  extrinsic  silicon  require 
cooling  to  liquid  helium  temperatures. 


Infrared  laser  frequency  mixing,  and  tunnel  time  measurements  using  laser 
rectification  in  SIM  junctions  have  offered  the  opportunity  to  explore 
variable  gap  metal  on  metal  and  metal  on  silicon  tunnel  devices  in  a 
controlled  way.  The  sharp  or  pointed  tip  geometry  introduces  a  geometric 
asymmetry  into  the  tunnel  junction  which  produces  a  field  enhancement 
for  tunneling  from  the  emitter  tip.  This  tunneling  asymmetry  provides  a 
mechanism  for  rectifying  an  alternating  laser  field. 

GANT  Fabrication 


Atomically  sharpened  silicon  emitters  are  formed  from  a  heavily  doped  n- 
type  substrate.  The  emitters  form  small  field  emission-like  structures  that 
are  positioned  within  30  Angstroms  from  a  planar  substrate  by  growing  a 
thin  tunnel  oxide  on  the  silicon  emitter.  A  thicker  deposited  oxide  is  used  to 
planarize  the  surrounding  emitter  surface.  One  of  the  drawbacks  to 
fabricating  nano-sized  junctions  with  macroscopic  integrated  circuit 
processes  is  the  device  variability  which  is  introduced  as  the  feahure  size 
becomes  progressively  smaller.  Wet  chemical  etch  as  well  as  reactive  ion 
etching  of  the  thick  planarizing  oxide  opens  windows  in  the  oxide  which 
are  subject  to  variation  due  to  the  diffusive  nature  of  the  etching  process. 
The  opening  in  the  thick  oxide  effectively  determines  the  emission  area  for 
each  of  the  field  emitters  in  the  array.  The  critical  spacing  between  anode 
and  cathode  is  unaffected  since  that  critical  distance  is  set  by  the 
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subsequent  growth  of  the  tunnel  oxide  atop  the  exposed  emitter  tip.  The 
fabrication  and  the  electrical  testing  of  the  devices  is  presented  in 
Appendix  A 

A  statistical  analysis  of  the  emission  area  variability  is  important  in 
determining  the  total  emission  current  developed  by  the  tunnel  array.  This 
realization  led  to  an  investigation  of  the  three  important  factors  which  will 
affect  the  final  design  rules  for  the  fabrication  of  GANT  devices. 

The  factors  are: 

1)  The  three  dimensional  nature  of  the  emitter  field  in  the  tunnel  junction. 

2)  The  effective  emission  area  which  may  be  significantly  different  from  the 
oxide  hole  opening. 

3)  The  statistical  variability  of  the  tip  radius  formed  by  stress  inhibited 
thermal  oxidation. 


Three  dimensional  tunnel  field 


In  all  previous  calculations  of  the  electric  field  in  an  STM  junction  we  have 
used  the  prolate  spheroidal  model  and  calculated  the  on-axis  electric  field. 
This  effectively  reduced  the  calculation  of  the  electric  field  in  the  tunnel 
junction  to  a  one-dimensional  case.  At  the  time  this  was  appropriate  since 
the  geometrically  correct  on  axis  solution  was  obtained  analytically  it 
represented  the  maximum  tunnel  junction  field  and  was  subsequently  used 
in  the  WKB  approximation  to  calculate  the  maximum  tunnel  current 
density. 

The  analytic  solution  for  the  full  three  dimensional  field  configuration 
including  a  fully  three  dimensional  image  and  multiple  image  interactions 
had  remained  an  unsolved  problem.  Although  numerical  methods  have 
been  suggested  these  solutions  suffer  from  the  usual  problem  of  accuracy 
without  a  comparison  to  an  analytic  solution.  In  Appendix  B  the  details  and 
results  of  that  calculation  are  presented. 

The  importance  of  this  solution  for  tunneling  in  geometrically  asymmetric 
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junctions  should  not  be  underestimated.  With  a  knowledge  of  the  spatial 
dependence  of  the  electric  field  reasonably  accurate  estimates  of  emission 
areas  are  possible.  These  results  are  discussed  in  Appendix  B. 


Statistical  Analysis  of  Emitter  Tip  Variation  and  Effect  on 
Current 


The  formation  of  atomically  sharp  and  nearly  atomically  sharp  field 
emitters  in  silicon  rests  on  a  technique  known  as  stress  inhibited  oxidation. 
This  oxidation  process  and  well  as  all  other  oxidation  processes  rely  on  the 
diffusion  of  the  oxidizing  species.  For  nearly  atomically  sharp  tips  the 
critical  oxidation  area  is  confined  to  a  region  less  than  100  X  100  Angstroms 
in  size.  Fluctuations  in  the  oxidizing  agent  concentration  produce 
variations  in  the  final  emitter  shape.  On  the  average  this  variation  may  at 
first  seem  negligibly  small.  In  fact  this  variation  has  a  exponentially  large 
effect  of  the  final  current  density. 

Unlike  the  case  of  planar  tunnel  structures  and  isoplanar  integrated 
circuits  these  devices  follow  a  Fowler  Nordheim-like  current  voltage 
characteristic.  The  exponential  nature  of  the  current  density  on  the  tunnel 
field  is  reflected  directly  in  the  sharpness  of  the  respective  tip  radii.  A  small 
variation  in  the  tip  radius  produces  a  large  change  in  the  emitter  current. 

In  the  analysis  we  have  performed  we  have  analytically  derived 
probability  distribution  functions  for  arrays  of  emitter  devices  for  two 
cases:  any  Gaussian  distribution  of  tip  radii  and  any  Rayleigh  distribution 
of  tip  radii.  The  controlling  current  voltage  relationship  was  assumed  to 
follow  a  classical  Fowler-Nordheim  law.  This  result  is  the  first  time  that  an 
analytic  probability  distribution  function  has  been  obtained  for  arrays  of 
field  emitter  devices.  This  result  becomes  an  important  analytic  tool  in 
characterizing  process  variation  and  predicting  array  performance. 


These  results  are  presented  in  detail  in  Appendix  C  and  are  generally 
applicable  to  field  emitter  arrays  independent  of  the  tip-anode  spacing.. 
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Inelastic  Processes  in  Tunneling 


The  shape  of  the  tunnel  current  as  a  function  of  tip-anode  spacing  is 
determined  using  kinetic  theory  and  the  transfer  Hamiltonian  method.  The 
photo-assisted  inelastic  tunneling  current  is  given  in  terms  of  laser 
frecjuency,  electric  field  ( which  is  calculatable  from  the  exact  three 
dimensional  solution  given  in  appendix  B),  and  appropriate  barrier 
parameters  for  a  materially  asymmetric  junction.  The  calculation  of  the 
laser  stimulated  inelastic  current  is  based  on  a  free  electron  model.  The 
technique,  however,  is  extendable  to  include  the  effects  of  band  structure. 

The  results  of  this  calculation  are  presented  in  Appendix  D. 


Digitally  Controlled  STM _ _ _ 

A  detailed  description  of  the  computer  controlled  STM  is  presented  in 
Appendix  E.  The  STM  itself  is  directly  operated  from  remote  sites  through 
the  imiversity  network.  Remote  operation  proved  to  be  necessary  given  the 
location  of  the  Nanostructures  and  Materials  Laboratory  on  the  7th  Floor  of 
the  Engineering  Building.  The  location  of  the  laboratory  on  the  7th  floor  is 
adjacent  to  the  building  elevator  shaft  which  proved  a  regular  source  of 
vibrations  which  were  damped.  Additionally  the  build  is  adjacent  to  a 
trolley  line  which  produces  significant  low  frequency  noise.  Successful 
operation  of  the  STM  and  the  need  to  accurately  control  the  tip-anode 
spacing  required  that  virtuaUy  all  experimentation  be  performed  between 
midmght  and  six  am  in  the  morning.  The  experiments  are  set  up  during  the 
day  and  the  data  acquisition  can  be  accomplished  from  adjacent  offices  and 
via  long  distance  telephone.  Without  the  remote  accessing  of  the  system  the 
noise  characteristics  of  the  building  would  have  dominated. 

A  discussion  of  the  computer  system  designed  and  used  through  out  the 
grant  t  appear  in  Appendix  E. 
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Surface  Modification 


The  determination  of  the  operational  tunnel  time  in  a  laser  irradiated  STM 
requires  a  controlled  tip-base  spacing.  During  the  experimental  calibration 
of  the  system  in  preparation  for  the  laser  experiments  an  unexpected 
modification  of  the  base  surface  was  detected.  The  gold  base  became  pitted 
giving  rise  to  spurious  tip-base  spacings.  This  modification  of  the  surface 
did  not  occur  when  the  STM  was  operating  in  the  imaging  or  normal 
turmeling  mode.  This  mode  is  defined  as  tip-base  spacings  such  that  the 
STM  rV  characteristic  is  linear-  that  is  spacings  less  than  10  Angstroms. 
Surface  modification  occurred  only  when  the  STM  operated  in  the  Fowler- 
Nordheim  region-highly  non-linear  current  voltage  characteristic,  and  only 
for  the  case  when  the  tip  was  positively  biased.  This  condition  corresponds 
to  the  forward  bias  case  during  laser  irradiation. 

Preliminary  results  indicate  the  creation  of  circular  and  elliptical  pits  is  not 
related  to  heating  of  the  gold  anode.  At  the  time  of  this  report  preliminary 
results  indicate  that  the  creation  of  the  nano-size  pits  is  a  field  effect. 
Additional  temperature  calculations  and  experiments  to  determine  the 
mechanism  for  the  pit  creation  are  currently  underway  in  an  effort  to 
ascertain  the  formation  mechanism. 

These  results  on  surface  modification  during  static  IV  testing  of  the  junction 
came  about  due  to  an  upgrade  we  developed  for  our  imaging  software. 
Since  the  detailed  structure  of  the  surface  is  important  in  maintaining 
constant  tip-base  spacings  planar  anode  regions  are  sought  for  IV  testing. 
Once  a  region  is  selected  it  is  extensively  mapped.  Static  current  voltage 
measurements  are  performed  with  re-imaging  before  and  after  each  FV  is 
taken.  This  is  a  tedious  and  time  consuming  procedure.  During  such  testing 
the  field  induced  pits  were  discovered. 

An  important  series  of  experiments  are  about  to  be  started  which  require 
further  modifications  to  the  software.  Since  this  surface  modification  is 
believed  to  be  a  field  induced  effect-  the  transient  or  induced  current  may 
play  the  dominant  role  in  creating  the  pits.  Induced  currents  occurring 
when  the  system  is  switching  voltages  may  be  further  reduced  by  changing 
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the  number  of  incremental  voltage  steps  and  hence  eliminate  the  pitting. 

On  the  other  side  of  the  issue  is  the  important  observation  that  surface 
modification  in  the  STM  junction  is  related  to  field  and  not  temperature 
effects  could  have  a  significant  impact  of  nano-scale  surface  modification. 

Figures  of  nanoscale  field  induced  pits  are  shown  in  Appendix  F. 

Laser  Induced  Thermoelectric  Effects  in  STM  Junctions 

We  have  calculated  the  spatial  and  temporal  temperature  distributions  for  a 
conical  metal  emitter  under  laser  irradiance.  The  full  three  dimensional 
heat  diffusion  equation  is  solved  simultaneously  wit  h  the  Fourier  equation 
for  heat  flux,  assuming  no  radiation  losses.  Both  Joule  heating  and 
Thomson  thermoelectric  effects  were  studies. 

The  details  of  the  Green's  function  method  used  in  the  calculation  are  given 
in  Appendix  G. 
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Appendix  A 

Fabrication  and  Testing  of  Geometrically  Asymmetric 
Tunnel  Structure 


Geometrically  Asymmetric  Nano-timnel  Structures 
D.C.Schweitzer,  T.E.  Sullivan 


Department  of  Electrical  Engineering 
Temple  University,  Philadelphia,  PA  19122 
R.B.  Marcus 


Department  of  Chemistry 

New  Jersey  Institute  of  Technology 

Newark,  N.J. 

Paul  H.  Cutler 
Department  of  Physics 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


Introduction 

In  r^ent  years  there  has  been  renewed  interest  in  the  use  of  vacuum  microelectronic 
devices  for  flat  panel  displays  and  rf  generation.  In  this  paper  we  extend  those 
application  to  the  development  of  high  speed  optical  and  infrared  detectors.  The 
development  of  bo  A  the  point  contact  diode  and  the  STM  tunnel  junction  have  shown 
that  novel  geometrically  asymmetric  nanotimnel  device  structures  can  be  used  as  an 
alternative  method  of  detection  for  high  speed  optically  modulated  signals.  One  of  the 
most  promising  device  structiures  under  investigation  are  the  at  atomically  sharpened 
silicon  field  emitters.  While  these  devices  are  similar  to  the  electrochemically  etched 
metal  emitters  used  in  both  the  scaiming  tunneling  microscope  and  field  ion 

microscopes,  atomically  sharpened  silicon  emitters  have  many  advantages  over  these 
devices.  ° 

We  present  the  results  on  the  development  of  a  novel  asymmetric  tunnel  structure 
whi^  can  be  directly  incorporated  into  a  silicon  integrated  circuit  rasing  standard 
semiconductor  processing  technologies.  This  device  structure  is  referred  to  as  a 
Geometrically  Asjmunetric  Nanotunnel  Structure  or  GANT  structure.  The  GANT 
structure  consists  of  an  atomically  sharpened  silicon  field  emitter  separated  from  a 
planar  metal  electrode  by  a  thin  ( d<  30  Angstroms)  tunnel  gap.  It  will  be  shown  that  an 
emitter  processing  scheme  such  as  this  facilitates:(i)  the  fabrication  of  large  arrays  of 
identical  field  emitter  elements,  (ii)  the  fabrication  of  controllable,  repeatable  cathode  to 
^ode  contact  geometries,  and  (iii)  the  incorporation  of  field  emission  devices  into 
integrated  circuits  using  current  semiconductor  processing  technologies. 

In  addition,  electrical  (I-V)  characteristics  show  a  current  density  enhancement  of 
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vacuum  tunnel  gap,  into  an  integrated  circuit  has  been  greatly  simplified  through  the 
use  of  a  Aermally  ^own  silicon-dioxide  barrier.  The  use  of  a  low  temperature,  thermally 
grown,  silicon  dioxide  film  has  been  shown  to  provide  an  accurate  and  repeatable 
method  for  forming  thin  tunnel  gap  spacings  <  50  Angstroms. 

Although  the  GANT  structure  presented  in  this  paper  differs  in  many  ways  from  an  STM 
tunnel  junction,  it  has  been  shown  that  this  design  scheme  has  a  number  of  advantages 
over  STM/point  contact  hmnel  structures.  These  advantages  include  (i)  fabrication  in 
large  arrays,  (ii)  fabrication  with  a  known,  controllable,  and  repeatable  contact  geometry 
and  (iii)  can  be  readily  incorporated  into  standard  integrated  circuits  using  current 
processing  technologies. 


Tunnel  Oxide 
(31  A) 


CVD  Oxide 
(3000  A) 


Field  Oxide 
(1150  A) 


Metalizatlon 
(200  A) 


n+  Si 


n+  S! 
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Three-dimensional  electrostatic  potential,  and  potential-energy  barrier, 
near  a  tip-base  junction 

Lo-Hong  Pan  and  Thomas  E.  Sullivan 

Department  of  Electrical  Engineering,  Temple  University,  Philadelphia,  Pennsylvania  19122 

Vallorie  J.  Peridier 

Department  of  Mechanical  Engineering,  Temple  University,  Philadelphia,  Pennsylvania  19122 

Paul  H.  Cutler  and  Nicholas  M.  Miskovsky 

104  Davey  Lab,  The  Pennsylvania  State  University,  University  Park,  Pennsylvania  16802 
(Received  28  March  1994;  accepted  for  publication  17  August  1994) 

The  geometry  of  an  atomically  sharp  or  nearly  atomically  sharp  tip  in  proximity  to  a  planar  anode 
may  be  closely  approximated  in  the  prolate-spheroidal  coordinate  system.  An  exact 
three-dimensional  electrostatic-potential  solution  for  a  free  charge  in  such  a  tip/base  junction  is 
given  in  this  letter,  including  calculations  for  both  the  symmetrical  on-axis  case  a^  the  asymmetric 
off-axis  case.  An  exact  solution  for  the  potential-energy  barrier  b  also  given;  this  solution  has 
immediate  applications  in  three-dimensional  tunneling  studies  and  in  calculations  of  el«;tron 
trajectories  in  micron-  and  submicron-sized  field-emitter  arrays.  G  1994  American  Institute  of 
Physics. 


We  report  that  an  exact  three-dimensional  solution  to  the 
Laplace/Poisson  equation  has  been  obtained  in  the  prolate- 
spheroidal  coordinate  system.  This  result  b  significant  be¬ 
cause  the  prolate-^beroidal  coordinate  system  provides  an 
accurate  geometrical  ie{X'esentation  of  an  emission  tip  close 
to  a  planar  anode,  a  geometry  which  arises  in  scanning  tun¬ 
neling  microscope  devices  (STM)  and  related  tuimel  struc- 
turs.  In  tbb  letter  the  first  off-axis  calculation  for  the  poten¬ 
tial  due  to  a  free  electron  in  such  a  voltage-biased  tip/base 
junction  is  givetL 

Theoretical  understanding  of  quantum-electronic  phe¬ 
nomena  In  nanoscale  devices  b  critical  in  the  continued  de¬ 
velopment  of  nanoelectronics.  As  important  class  of 
quantum-electrooic  devices  consist  of  a  biased  tip/aiKide  or 
tip/base  junction,  in  which  the  tip  b  atomically  sharp  ir^<2 
nm)  and  dose  (3  run)  to  the  planar  surface.  For  example, 
such  devices  UKhide  scanning  turmeling  microscopes  (STM) 
and  geometrically  asymmetric  nanotunnel  structures 
(GANTS).*  Of  particular  interest  b  the  capability  to  quantify 
the  tumbling  current  for  such  devices,  a  calculation  which, 
in  turn,  requires  an  accurate  estimate  of  the  induced  potential 
energy  of  a  point  diarge  in  the  tip/base  fraction.  Much  of  the 
theoretical  work  to  date  has  been  carried  out  using  an  effec¬ 
tively  one-dimeostonal  or  planar  analysis,  and  whQe  there 
are  some  devices  (e.g.,  quantnm-weO  beterostructures)  for 
which  thb  one-dimensional  ^iproach  proves  to  be  an  ad¬ 
equate  strategy,  h  b  dearly  not  an  ideal  geometrical  approxi¬ 
mation  for  a  three-dimensional  device  utiliring  a  tip/base 
junction.  However,  h  emerges  that  the  geometrical  features 
of  devices  such  u  the  STM  may  be  closely  approximated  in 
the  pfolate-spherolda]  coordinate  system,  with  the  tip  and 
base  each  represented  as  a  hyperboloid  surface. 


energy  barrier.  It  b  worthwhile  to  note  that  exact  solutions 
such  as  these  permit  the  calculation  of  these  fields  to  arbi¬ 
trary  accuracy,  even  in  regions  of  severe  spatial  variation; 
calculations  are  given  here  f<»  the  on-  and  off-axb  case  of  a 
point  charge  in  a  biased  tip/anode  junction. 

To  begin  the  analysb,  the  prolate-spheroidal  coordinate 
system  b  given  by^ 

x=fl  sinh(tf)sin(i;)cos(^), 

y=a  sinh(u)sin(i;)sin(<^),  (1) 

z=a  cosh(u)cos(t;), 

where  a  b  one-half  the  dbtaiKe  between  the  hyperboli  fod; 
curves  corresponding  to  constant  v  are  hyperboloid  surfaces, 
while  those  for  constant  u  are  prolate-s^roidal  surfaces. 
Geometrically  the  problem  has  axial  symmetry,  and  if  the 
charge  b  on  the  symmetry  axb  «»0  the  solution  for  the 
resulting  electrostttic  potential  b  effectively  a  one¬ 
dimensional  problem  which  could  be  addressed,  say,  with  a 
method-of-images  technique.’-*  However,  if  the  change  b  lo¬ 
cated  at  an  off-axb  position  then  the  problem  for  the  electro¬ 
static  potential  b  fully  three  dimensional. 

Starting  with  (1),  the  variable  transformations 
{=eosh(a),  i7=cos(o),  generate  a  coordinate  system  defined 

by 

1  Vl  -  ^  C08(^), 

y=a>/f^-l  Vl-i?’  sin(^),  (2) 

The  region  between  the  two  hyperboloid  surfaces  Is  the 


(_1)»  r(i+iT+in) 

TtMh(irT)  r(|+iT-m) 


X^t-t').  (10) 

note  that  the  second  sotedon  for  is  ftom 

(8)  it  is  evident  that  the  solutions  f«-  5(17)  are  also  conical 
functions.  Consequently,  the  general  solution  for  V|,  noting 
that  the  eigenvalues  m  are  discrete  and  the  eigenvalues  rare 
continuous,  is 
«• 

V;=2  (2-«o.«)cos(m^). 

■1-0 

f- 

I  t*7*?(v)+CKT(->»)l*T(£)<<T. 


with 


k7(  -  vc)K(  vt)  -  vi) 

Here,  the  charge  is  located  at  Xe»({,.^,,0),  and  the 
constant  a?  in  (12)  is 


^  rtanh(irT) 
cosh(irT) 


r(i+fr-»n)  * 
r(J+iT+m) 


(13) 


Thus  the  electrostotic  potential  V  (4)  is  calculated  from 
(6),  (7),  and  (11),  and  the  potential  energy  d>  is 


b:=—  a:K:u,)K':{- m)- 

4eo« 

7<)C(  -  7i)  -^7(  -  7i ) 
K:{vx)K7(-m)-f(ri-  >7i)^T(»72)  ’ 


4>=9(V,+  V;), 


(14) 


(12) 


where  Vg  is  given  by  (6)  and  V,  by  (11).  Also,  it  can  be 
shown’  that  the  electtic  field  on  the  emitter  surface  at  rj-  77, , 
El, is 


-2Vo 


\  1-J72  1  +  7,/ 


•  X  (-l)"(2-^o.J 


xcos(m«)|^r|  ,unb(„T)  1JC({,)*T(£)]|  1 


r(i+iT— m) 


k:{  -  7c)/C(  72)  -  7c)i^7(  -  72)  \ 


r(i+iT+in) 


(15) 


Note  that  the  first  tenn  in  Eq.  (15)  is  the  eminer-surface 
electric  field  due  to  the  bias  potential  which  is  independent 
of  the  charge  location;  the  second  term  is  the  interaction  field 
produced  by  the  induced  surface<harge  distributions  on  the 
tip  arid  base  and  the  frce-charge  located  at  7,,0). 

Figure  2  shows  computed  contour  plots  of  the  potential 
V  due  to  a  point  charge  in  a  voltage-biased  tip-base  junction, 
for  both  the  on-axis  and  an  off-axis  case.  In  both  plots: 
rji  =0.9  (which  correqioods  to  a  tip  half-angle  of  0.45  radi¬ 
ans  or  26*);  the  value  of  o  is  chosen  so  that  the  tadhis  of  the 
tip  at  the  apex  is  20  A;  and  the  bias  voltage  of  the  tip  relative 
to  the  base  is  5  V. 

In  conclusion,  an  exact  three-dimensional  solution  for 
the  electrostatic  potential  doe  to  a  charge  in  a  voltage-biased 
tip/anode  junction,  modeled  in  the  prolate-spheroidal  comdi- 
nate  system,  has  been  obtained.  The  results  of  this  study  will 
be  used  for  characterizing  the  tunneling  current  at  a  tip^anode 


junction  using  a  three-dimensional  tunneling  analysis,  and 
will  be  used  to  cany  out  a  systematic  study  of  off-axis  tun- 
and  its  impact  on  the  estimation  of  active  emission 
areas  in  nanosized  devices. 
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(k) 
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•®^  •'*  •  **  poteatUJ  energy  due  to  the  tip/anode  *»»««  (e)  Th 
rubtidUry  p«,ble«  far  .  fl*  poteSel  en*,^  doeVdSd*^. 


VJ*. 


«*(!*- 7^) 

e-r? 


(f’-l)(l-7*)^**J 


(3) 


The  Laplacian  operator  is  linear,  and  consequently,  the  solu¬ 
tion  for  the  potential  V  can  be  expressed  as  the  sum  of  two 
subsidiary  potentials, 


V=V,+V, 


(4) 


Here,  V,  is  the  potential  due  to  the  bias  voltt^  applied  to 
the  tip-plate  junction.  f<x  instance  Vq,  and  includes  the 
effects  of  the  chai:ge’s  interaction  with  the  sutface- 

chaije  distribttdoo  on  the  base  and  tip;  note  that  this  surface 
change  is  distributed  so  as  to  yield  zero  potential  on  these 
conducting  surfaces.  Schematics  associated  with  V  and  the 
subsidiary  problems  V, ,  are  given  in  Fig.  1. 

From  Fig.  1  (b)  it  is  evident  that  the  problem  for  Vg ,  the 
potential  due  »lely  to  the  potential-bias  between  the  tip 
and  aiKxie,  is  independent  of  f,  and  the  i^ysical  location 
of  the  charge  x, .  Thus,  the  equations  governing  Vg  are 


<?  f  ,  d 


Vg-O,  i'a(7i)  =  i'ot  i'a(72)=0. 


(5) 

solution  for  Vg  may  be  found  by  elementary  mr-jhnfjc 
-nd  is 


l(^)(S)l 


(6) 


ie  potential  doe  to  the  diarge  satisfies  the  Pbissoneqoa- 
n  and  may  be  written  as  the  sum  of  a  general  (homoge- 
ous)  arxi  a  partkolar  solution,  viz. 


4ir€o|*-x,| (7) 

^  interaction  potential  V,  satisfies  the  Laplace  equation 
has  homogeneous  boundary  conditions;  however;  shme 
potential  it  fuUy  three  dimensional  it  is  cooveniem  to 
karate-  the  solution  into  a  product  of  functioos  in  each 
dmatc  direction,  for  instance  »',=«(f)5(i7)F(^).  Ap- 
ig  the  Laplacian  operator  (3)  on  K,  in  this  separated 
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form  and  designating  the  separation  constanu  as  A  and  m^, 
the  solutions  in  each  coordinate  variable  are  governed  by 


To  obuin  F{4>)  the  charge  may  be  set,  with  no  loss  of 
generality,  to  an  azimuthal  location  ^0,  so  that  the  solu¬ 
tions  for  F(^)  will  be  of  the  form  cos(m^). 

TV  boundary  conditions  for  V,  are  homogeneous  at  i* , 
ffi't  ihis  suggests  that  R{Q  must  be  expanded  in  an  orthogo¬ 
nal  ^  of  functions  and  thus  the  constant  A  must  be  negative. 
S^g  A =n(n  + 1)  gives  associated  Legendre  polynomials 
P,  with  complex  degree  n*-j±/r.  Following  Hobson.’ 
designate 

the  functions  K"  are  called  conical  fiitxrtions,  and  the  or¬ 
thogonality  relation  for  the  conical  functions  is* 


Mh  ftl  “******  rymaeuy 

«i^WCoiiio«  of  eoQftamiKMwtitJ 


(_1)"  r(i+iT+m) 

tanh(irT)  r(^+iT-m) 


X5(r-r').  (10) 

note  that  the  second  sobtion  for  K^iji)  is  ftom 

(8)  it  is  evident  that  the  solutions  for  5(17)  are  also  conical 
bncdons.  Consequently,  the  general  solution  for  V; ,  noting 
that  the  eigenvalues  m  are  discrete  and  the  eigenvabes  r  are 
continuous,  is 


ai-0 


+  V)]f^riOdT, 


with 


4«o« 

’72)-0-  Vi)f(7iV2) ' 


(11) 


(12) 


C:=^aWCc)^(i7,). 

^:(-i7c)C(’72)-^7(>7c)^^7(-  172) 

K:ivi)K:{-v2)-K';i-vi)K:{v2)  ■ 

Here,  the  charge  is  located  at  x,=(^(,i7j,0),  and  the 
constant  <*7  in  (12)  is 


a7= 


T  tanh(irT) 


cosh(irT) 


r(i+ir-m) 


r(|+iT+m) 


(13) 


Thus  the  electrostatic  potential  V  (4)  is  calculated  from 
(6),  (7),  and  (11),  and  the  potential  energy  <I>  is 


<*»=q(V,+V;), 


(14) 


where  Vg  is  given  by  (6)  and  V,  by  (11).  Also,  it  can  be 
shown’  that  the  electric  field  on  the  emitter  surface  at  , 
£1,  is 


-2V, 


1 


£,=■ 


-  .  ■  — - + - p==^7=-S  (-1)"(2-5«.J 

"  >/f’-»7?Vl->7?i  +  2irflfo  Vf’-7iVl-’7?  "-0 


1-J72  1  + J7i) 
r(^+iT— m) 


Xcos(m<^)  I  dT  Ttanh(rrT) - [£7(fe)£^7(f)] 

r(^+i'T+m) 


£7(  -  Vcm  >72)  -£:(  ’7c)£T(  -  V2) 
£7(  ni)K':(  -  V2)  -£:(  -  >7i  )£■:(  »72) 


(15) 


Note  that  the  first  term  in  Eq.  (15)  is  the  emitter-surface 
electric  field  due  to  the  bias  potential  which  is  independent 
of  the  charge  location;  the  second  term  is  the  interaction  field 
produced  ^  the  induced  surface-charge  distributions  on  the 
tip  and  base  and  the  fiee-charge  located  at  Xe=(fc»’7e*0)- 

Figure  2  shows  computed  contour  plots  of  the  potential 
V  due  to  a  point  charge  in  a  voltage-biased  tip-base  junction, 
for  both  tte  on-axis  and  an  off-axis  case.  In  both  plots: 
rjy  =0.9  (which  corresponds  to  a  tip  half-angle  of  0.45  radi¬ 
ans  or  26**);  the  value  of  a  is  chosen  so  that  the  radius  of  the 
tip  at  the  apex  is  20  A;  and  the  bias  voltage  of  the  tip  relative 
to  the  base  is  5  V. 

In  conclusion,  an  exact  three-dimensional  solution  for 
the  electrostatic  potential  due  to  a  charge  in  a  voltage-biased 
tip/anode  junction,  modeled  in  the  prolate-spheroidal  coordi¬ 
nate  system,  has  been  obtained.  The  results  of  this  study  will 
be  us^  for  characterizing  the  tunneling  current  at  a  tip/anode 


junction  using  a  three-dimensional  tumieling  analysis,  and 
will  be  used  to  carry  out  a  systematic  study  of  off-axis  tun¬ 
neling  and  its  impact  on  the  estimation  of  active  emission 
areas  in  nanosized  devices. 
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Abstract 

Classical  field  emission  theory  is  based  on  the  one-dimaisional  planar  Fowlo’-Nordhdm 
theory  which  relates  current  density  to  die  applied  constant  electric  field  at  die  onitter 
surface.  However,  in  the  actual  non-planar,  dp-anode  geometry  die  dectric  fidd  varies 
along  die  microtip  surface,  giving  rise  to  a  spatially  dependent  currmt  density.  This 
prqier  accounts  for  the  electric  field  and  die  current  density  variations  aloi^  a  miootip 
and  derives  analytical  expressions  fOT  tip  current  and  emission  area.  CalculationB  for 
current  density  and  total  tip  current  are  carried  out  using  an  exact  diree>dimaisiQnal 
dectric  field  solution  in  prolate  spheroidal  coordinates  in  coryunction  widi  die  dassical 
FowlCT-Nordheim  emission  dietny.  Die  numerical  results  provide  a  pttysical  context  for 
dM  oftoi  used  fidd  enhancement  fiu;tor  p  through  die  use  of  exact  analytical  expressions. 
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requires  an  exact  tfuee-dimensional  sdution  for  tfie  dectrostatic  potailial  both  at  die 
emitter  surface  and  in  die  adjoining  free  space  region  between  emitter  and  anode. 

One  approach  often  taken  for  ad^ting  the  one-dimensional  Fowler-Nordheim  analysis 
to  inherendy  three-dimensional  structures,  is  to  model  die  surface  electric  field  as  some 
constant  times  die  applied  voltage  V,  i.e.  F=pV.  Here,  the  parameter  p  or  field 
enhancement  factor  is  often  modelled  as  a  function  of  the  tip  geometry.  For  exanqile, 
one  model  for  the  field  enhancement  &ctor  is  P=l/kr,  where  r  is  the  microtip  radius  of 
curvature  and  k  is  a  factor  to  be  empirically  determined  [2].  p  factors  are  routinely 
reported  which  provide  reasonable  fits  of  cunent/Voltage  data  to  die  Fowler-Nordhdm 
equation,  assuming  F=pV  [2,3].  Alternative  apfmiaches  have  also  been  tried  to  fit 
Spindt-emitter  data  using  a  model  which  involved  two  empirically  determined  factrxs, 
p  and  y,  to  relate  die  electric  field  quantity  F  to  die  gate  and  collector  voltages, 
respectively  [4]. 

Theoretical  difidculties  for  the  field  enhancement  fKtor  are  known  to  arise  in  the  context 
of  extiemely  sharp  tips  (r<100  A).  For  exanqile,  die  Fowler-Nordheim  anafysis  has 
built-in  assumi^ons  concerning  die  electron  8iq>[dy  function  and  induced  surface  charge, 
assunqidons  which  are  predicated  on  a  planar,  semi-infinite  medium  emission  sur&ce; 
for  die  case  of  extremely  sharp  tips  eidier  die  radius  of  curvature  or  die  physical 
dimension  of  the  tip  may  violate  such  assunqitions  [S].  An  additional  conc^itnal 
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difficulty  with  tfie  field  enhancement  q>proadi  is  fiiat  goienlly  file  field  is  treated  as 
a  constant  over  the  entire  emissioa  surface.  Recognizing  fiiis,  a  more  recent  fiieoretical 
analysis  has  endeavored  to  quantify  p  by  modeling  the  q)ex  of  the  microtip  as  a  floating 
sphere  [6]. 

A  principal  restriction  of  the  field  enhancement  concq>t  for  ad^’ng  die  Fowler- 
Nordheim  to  diree  dimensional  structures  lemaiiis;  that  is,  a  calculation  of  tuimel  current 
cannot  be  obtained  a  priori  for  a  given  {rfQ^sical  system  because  an  oiqnrical  estimate 
of  p  is  required.  Thus,  alfiiough  the  p  futor  provides  a  cmivaiient  scheme  for  die 
comparison  of  experimental  data,  it  is  strict  an  enqurically  determined  quantify.  What 
is  desired  ben  is  a  conceptual  fiamewock  wMdi  permits  an  estimate  of  tunnel  current, 
in  the  context  of  die  Fowler-Nmdheim  anafysis,  fitxn  first  prindides;  diat  is  die 
objective  addressed  in  this  paper. 

The  work  described  in  dtis  pqier  offers  a  scheme  to  evaluate  the  performance  of  field 
onission  microtips  firom  first  principles,  in  the  contact  of  die  Fot^er-Nordheim  dieory, 
widiout  the  need  of  an  enqnricalty  detomiiied  p  factor.  This  objective  has  bear 
achieved  by  die  use  of  an  exact,  fully  diree-difflaisional  sdudmi  for  die  dectric  fidd 
in  a  microtip/base  junction  [7].  This  geometry  modds  the  nucrodp  and  anode  as 
hyperboloid  surfiices,  and  obtains  an  exact  sotudcm  for  sevnal  fuOy,  threenfimensicmal 
dectrostadc  quantities  due  to  a  flee  charge  in  a  tqi/base  junction,  in  die  prolate- 
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spheroidal  coordinate  system.  To  fully  appreciate  die  significance  of  this  model  it  is 
wordiwhile  to  consider  Figures  1  &  2.  Figure  1  shows:  (a)  a  tip  modeled  as  a  constant- 
radius  structure  (such  as  that  idealized  in  the  floating-sidiere  models  die  dotted  line,  and 
(b)  die  tip  modeled  as  a  hyperboloid  surface,  die  solid  line.  In  Figure  2  a  transmission 
electron  micrograph  (TEM)  of  a  fabricated  silicon  tip  is  reproduced  [8],  overlaid  by  die 
two  models  ^en  in  Figure  1.  It  is  self-evident  diat  die  Itypetboloid  surface  oflFers  an 
excellent  physical  approximation  for  the  geometrical  features  of  a  field  emission 
microtip,  and  such  agreement  motivated  the  utilization  of  die  exact  prolate  spheroidal 
solution  for  the  electric  field  to  quantify  die  emission  characteristics  of  microtips. 

The  methodology  adopted  in  dus  p^qier  is  as  fdlows.  Ibe  exact  solution  for  the 
electrostatic  field  of  die  microtip/base  junction  given  by  die  prolate  spheroidal  model 
is  substituted  into  die  Fowler-Nordheim  equation,  yidding  an  analytic  expression  diat 
describes  the  tunnel  current  densify  at  ev^  point  on  die  surface.  The  total  tip  current 
is  then  obtained  by  integration  of  diis  current  densify  expression  over  the  surface  of  the 
microtip.  The  current  densify  on  the  tip  surface  decays  r^dly  widi  the  distance  firom 
the  tip  apex.  Using  this  result,  we  have  been  able  to  devdop  an  estimation  sdieme  for 
the  effective  emission  area  of  the  tip. 
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n.  Electric  Field  Derivadoii 

Ibe  doivation  for  die  electee  field  on  a  miootip  surface,  in  a  biased  tip/base  junction 
as  modeled  in  a  prolatd-ipheroidal  geometry,  is  given  in  [7].  Ibe  prolate  sphermdal 
coordinate  system  shown  in  Rguic  3  is  defined  ly  die  coordinates  ii,  and  ^  [7]  whoe 
\  e  [l,aX  T|  €  ^  [0,271].  til  COTcsponds  to  die  bypeiboloidal  suifKe  of  die 

dp,  til  is  die  planar  anodc  suifiu^e,  lines  of  constant  %  define  prolate  spheroidal  suifiues, 
and  is  die  ayimntfial  angle  of  rotadon  about  die  symmetry  or  major  axis.  Odier  key 
parameters  of  interest  include  the  fod  distance  of  2a,  die  dp-anode  distance  d,  and  die 
«ip  half  angle  ft  These  three  parameters  are  idated  by  dw  expression  e=co8*\d/a).  The 
dp  radius  r  on  axis  is  defined  as  r=a  smO  tan6. 

a  voltage  Vj  is  qiplied  betweai  die  tqi  and  anode,  the  Li^lace  equadem,  V*V=0,  must 
be  satisfied.  In  prolate  sidieroidal  coordnates,  die  Lrqilacian  is  givoi  by: 


Setting  V(iii)=0  at  the  dp  and  V(ii])=V,  at  the  anode,  die  dectrostadc  potential  is 
symmetric  in  %  and  <{>  and  dms  is  indqiendent  of  %  and  ^  so  that  die  ^udon  for  V  is: 
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(2) 


The  field  F(4)  at  any  point  on  die  surface  of  the  tip  is  determined  fi-om  the  negative 
gradient  of  V(ti).  F(^)  is  a  vectw  widi  components  in  the  and  (Ji  directions  equal  to 
zero.  The  resulting  magnitude  of  F(J^)  in  die  q  direction  is: 


am  (9)  VF^5?W  m(ii|||^) 


(3) 


where  cos(6)  has  been  substituted  for  q  [9]-  Note  that  6  is  the  cone  angle  of  die 
microdp,  widi  6->7i/2  a  very  flat  dp,  and  d->0  an  extremely  sharp  dp.  Note  also  that  \ 
in  effect  measures  die  distance  on  the  top  surface  fi'om  the  apex,  with  4=1  the  apex 
ordinate. 
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DL  Impiicatioai  for  Field-Enhancemeiit  Factor 

It  is  possible  to  substitute  an  expression  in  r,  tbc  radius  of  curvature  at  die  2q>ex,  for  die 
parameter  a  and  put  F(^)  in  a  form  that  follows  die  field  enhancement  concept  F=PV, 
where  p  is  expressed  as  p=l/kr. 


P(E)  -  ^0  -2  tan(e) 


(4) 


With  F(i|)  in  this  form  a  field  enhancement  factor  defined  by  k,  would  be  e:q>ressed  as: 


*<«■ - <=> 

Asymptotic  analysis  may  be  used  to  show  that  k»l  in  the  limit  of  very  sharp  tip  (6->>0), 
and  that  kr->d  in  the  limit  of  very  blunt  tips  (6->n/2,  ^-^1).  This  means  diat  as  die  tip 
Iffoadens  to  become  a  |danar  sutfiu:e;  the  geometry  appcan  as  two  paralld  {dates  ai»l 
the  p  fiu:tor  s^iproaches  the  ex{>ected  limiting  value  lAL 
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IV.  Tip  Current  Calculation 

The  most  general  form  of  die  Fowler-Nordheim  equadem  is  [10]: 


3  IT* 


8x2n|rC 


exp 


1  3J2eF  ' 

\  ♦  JJ 

(6) 


In  dus  expression  F  is  the  electric  field  (V/cmX  v{/  is  the  imUerial  work  fimetion  (eV), 
h  is  Planck's  constant  (og-seeX  e  is  the  electron  charge  (stat-coulX  j  is  die  current 
density  (anqj/cm^  and  m  is  die  electron  mass  (gm)> 

The  V  function  is  dependent  on  an  involved  competition  of  elliptic  int^rals  and  varies 
fiom  1  to  0  as  electric  field  strength  increases  [10].  The  t^  funedem  is  a  slowly  varying 
-function  of  die  applied  field  and  is  approximated  as  1.1  [10].  In  order  to  avoid 
calculations  with  the  int^rals  that  define  v,  an  approximation  that  has  been  made  is 
[11]: 


v{F)  ^  0.95-y{F)^  yiF)=2.79xlO'*^  (7) 
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Using  Ais  approximation  and  substituting  die  values  for  ^  and  die  constants,  die 
equation  can  be  expressed  in  CGS  units; 


(8) 


This  form  of  die  equation  uses  dectric  field  values  in  volt/cm.  Converting  die  field 
values  to  volt/A  the  equation  can  be  written  as: 

^(i^=AF*©xp|-^j  (9) 

where  the  constants  A  and  B  are  given  by: 

ji-  ^ •  54^0^° 9 .^1  j  S=0.649f'^  (10) 

The  Fowler-Nordheim  equation  develqied  above  provides  for  the  currmt  density  j(F) 
measured  in  amps/cmi^  fix’  die  case  of  a  (danar  geometry  and  a  constant  electric  fidd. 
In  comparison,  die  prolate  sidiermdal  dectric  field  F(^)  varies  widi  ^  die  dimoisionless 
distance  from  die  ^lex  of  the  nuCTOtip,  whoe  the  fidd  F  is  greatest 
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Unlike  die  one-dimensionai  analysis  where  die  field  is  c<a)stant  bodi  mi  die  surface  and 
in  the  adjoining  five  space,  die  three-dimaisional  nature  of  F(^)  yields  a  spatially 
variable  electric  field  bodi  on  die  surfece  and  in  the  adjoining  five  space.  If  we  imagine 
the  tip  to  be  conqiosed  of  numerous,  infinitesimal  planar  s^ments  of  area  d^d4> 
distributed  in  a  manner  similar  to  the  marks  on  a  golf  ball,  we  can  apply  the  Fowler- 
Nordheim  equation  on  each  planar  segment  and  then  integrate  die  resulting  j(^)  over  die 
microdp  sur&ce  to  obtain  tip  current,  where  4o  ^  ^  iqqier  limit  of  ^  diat  defines  die 
effective  emission  area;  die  systematic  selection  criterion  for  die  4o  cutoff  value  is 
described  below. 


j  (5) 


The  tip  current  i  is  exidicidy  shown  as  a  function  of  various  parameters,  and  the 
functions  h^  and  h^  are  prolate  spheroidal  metrics  [9].  Substitution  of  the  various 
functions  in  and  sin^lifying  leads  to  die  e?q;>ressi(m: 


ii  ^5*-cos*(0) 


(13) 
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where  the  functions  x  and  y  are  defined  as: 


a.sisexio**  e: 


x(8,t)  = 


if) 


(14) 


y(d,8,t)  =0.3245  dt^tan(e)lx>(ll|||<|j.] 

Numerical  results  based  on  die  above  method  have  been  selected  and  the  corresponding 
curves  for  electric  field,  current  density,  and  tip  current  are  plotted  in  Figures  i-6.  It 
is  aignifirant  to  nodce  how  quickty  the  current  density  is  reduced  as  a  function  of  ^  die 
distance  along  the  surface  fiom  the  tip  axis.  Various  sets  of  calculations  have  bera 
done  for  other  values  of  V,  d,  and  0,  mclutfing  die  limiting  gecxnetries  of  a  veiy  sharp 
tip  and  a  tip  so  blunt  that  it  approximates  a  flat  plane. 


In  those  cases  where  die  value  of  in  the  region  of  electrm  emission  remains  close  to 
1,  the  dp  current  integral  of  (13)  can  be  very  closely  approximated  as: 


i(d,V,0,f)=i^  [exp  (0) J-eocp  8ln(0)Jj  (15) 
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lliu  equaticm  provides  the  value  of  tip  current  in  amperes  for  a  ge(»netiy  specified  by 
d  and  6,  and  for  an  effective  emission  area  defined  by  die  range 


V.  EfTective  Emission  Area 

Calculations  for  the  tip  current  indicate  dial  most  of  die  current  is  due  to  die  tip  area 
that  is  very  close  to  die  tip  axis.  The  numerical  results  show  that  99.9%  of  the  tip 
current  is  fiom  an  area  where  current  density  has  decayed  to  3  orders  of  magnitude 
below  die  current  density  peak  located  on  axis.  This  lower  limit  of  3  orders  of 
magnitude  reduction  in  current  density  will  be  used  in  a  derivation  of  the  effective 
onission  area. 

The  value  of  ^  at  which  the  current  density  is  reduced  by  a  factm'  of  lO'^  is  given  by: 


(16) 


The  i«e-exponential  terms  for  j<4)  and  j(lX  fitm  equation  (IIX  dOfeedvdy  caned  for 
values  of  ^  close  to  1,  and  die  exponential  terms  determine  die  value  of  %  winch 
satisfies  (16).  Performii^  dns  (^>erati(»i  and  solving  for  ^  yidds  die  fdlowii^ 
expression: 
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-21.287  V 


•^sln(e) 


+C08*  (0) 


(17) 


Hds  expression  defines  file  proper  4o  ^  current  equation  shown  in  (IS).  In 

file  case  (b)  of  fiie  parameters  selected  for  fiie  figures,  the  value  of  ^  is  1.0068.  Wifii 
file  substitution  of  (17)  into  (IS)  and  simplifying,  the  equation  fin*  t^  current  can  be 
expressed  as: 


Similaiiy,  (17)  can  now  be  used  to  calculate  file  effective  emissimi  area: 

(19) 

A{d,e,  U)  -«d«8ln(e) 

-8  d*taii*(0)+  xd^8ln(d)  ln{l+8ln(6)] 

In  the  case  (b)  of  the  parameters  sdected  for  the  figures,  the  efifecfive  onission  area  is 
32,706  A*. 
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VL  Comparisons  wftfa  Experimental  Detenninations  of  Beta 

The  expression  for  a  local  P,  the  field  enhancement  fsctor  associated  widi  the  expression 
is: 


P(5)  = 


-2 


(20) 


The  local  beta  factor  is  dependent  only  upon  the  gemnetiy  of  die  tip  defined  by  the  tip 
half  angle  6  and  the  fod  distance  a,  and  varies  idong  the  surface  of  die  tip  widi  values 
of  Figure  7  exhibits  dus  behavirn*  by  displaying  die  inverse  of  die  beta  factor  versus 
for  a  range  of  sdected  on  axis  tip  radii.  Note  diat  die  value  of  P'*  converges  to  die 
parallel-fdate  value  of  the  tip-anode  distance  d  in  die  limit  as  the  tip  radius  approaches 
infinity. 

The  following  experimental  estimates  for  an  averaged,  global  tip  value  for  p,  as  given 
by  several  investigators,  are  in  close  agreement  widi  the  results  shown  in  Figure  7. 


Year 

Researcher 

Tip  Radius 

P  value 

Reference 

1953 

PykeTTrolan 

18  pm 

12,800  cm  * 

[12] 

1993 

Jenkins 

~10  nm 

636,000  cm  * 

[13] 

1994 

Hrmg  et  oL 

~25  nm 

210,000  cm  * 

[14] 
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Vn.  Conclusiom 


The  use  of  the  prolate  spheroidal  tip/anode  model  [7]  allows  for  the  exact  representation 
of  die  bias  electric  field  and  is  used,  in  coiyunction  with  classical  Fowler-Nordheim 
tiieoty  [1]  to  surface  varying  electric  field  tuimel  current  densities,  net  tip 

current,  and  effective  emission  area.  Because  the  prolate  spheroidal  coordinate  system 
is  such  a  close  geometrical  approximation  to  the  tip-to-plane  geometry,  our  microtip 
calculations  widun  that  system  possess  an  unexpected  sinq[>licity. 

It  has  been  suspected  for  years  that  the  p  factor  is  only  geometry  dependent,  and  data 
has  shown  it  to  be  of  constant  value.  The  exact  field  solution  in  prolate  spheroidal 
coordinates  largely  corroborates  those  assertions.  As  derived  for  a  single  nricrotip  die 
P  factor  is  solely  dependent  on  geometry;  it  does  not  vary  with  any  electrical  opoating 
characteristics. 


Close  examination  of  Figure  7  reveals  that  experimentation  with  microtips  of  radii  250- 
400  A  (e^.,  Spindt  tips)  would  observe  a  nearly  constant  p  fiwtor.  As  researchers  find 
ways  to  create  rruCTotip  arrays  widi  an  average  of  tip  radu  approadung  50  K,  dMy  may 
begin  to  observe  a  non-constant  value  for  p. 

In  die  derivation  of  the  prolate  sfdienndal  tip  currm^  die  exact  sdutirm  for  die  dectric 
field  F(^)  was  substituted  into  die  classic,  fdanar  FowlCT-Nordhdm  equation  for  current 
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density,  and  we  imagined  die  tip  as  numerous  planar  segments  of  area  dl^d(l>.  In  reality, 
as  die  dectron  tunnels  fixwn  the  surface,  it  does  not  experience  a  conqilctely  flat  surface 
as  assumed  by  the  Fowler-Nordheim  dieory.  However,  die  implicit  assumption  here  is 
that  the  surface  jqipears  approximately  planar  as  long  as  die  local  radius  of  curvature 
for  die  surface  is  larger  dian  die  tuimel  distance  of  the  electroiL 


17 


references 

'RJi  Good  Jr.  and  E.W.  Muller,  Handbuch  der  Phvsik,  XXI,  Si»inger,  192,  (1956). 

^rodie  and  C.A.  Spindt,  Vacuum  Microelectromcs^  Academic  Press  Inc.,  New  Yoik, 
14,  (1992). 

’WJ*.  Dyke,  Sd.  Am.,  210-1,  108,  (1964). 

^avid  W.  Jenkins,  TF.FF  Trans.  Electron  Devices,  ED-40,  666,  (1993). 

’Jun  He,  Theoiy  of  Electron  Emission  from  Atomically  Shaip  Metallic  Emitters  in  High 
Electric  Fields,  Uie  Pennsylvania  State  Univeraity,  hr  (1992). 

^an  Nicolaescu,  J.  Vac.  Sd.  Technol.  B  11,  392,  (1993). 

’LJl  Pan,  TE.  Sullivan,  VJ.  Peridier,  PIL  Cutler,  and  NAl  Miskovsky,  Appl.  Phys. 
Letters,  65  (17),  2151-2153  (1994). 

*PhotognqA  courtesy  of  RE.  Marcus. 

*G.  Aifken,  Mathematical  Mediods  for  Physicists,  Academic  Press  Inc^  New  Yoric,  103- 
106,  (1970). 

*®RJ1  Good  Jr.  and  E.W.  Muller,  188. 

”C.A.  Spindt,  I.  Brodfie,  L.  Humphrey,  and  EE.  Westeiberg,  J.  Appl.  Phys.  47,  5250, 
(1976). 

“WP.  Dyke  and  JJK.  Troian,  Physical  Review  89  (4X  805,  (1953). 

“David  W.  Jenkins,  670. 

Hong,  ML  Aslam,  M.  Feldman,  and  M.  Olinger,  J.  Vac.  Sa.  TechnoL  B  12,  767, 
(1994). 


18 


CAPTION  PAGE 


Figure  1.  Conq)arison  of  prolate  spheroidal  and  constant  radius  tip  geometries. 

Figure  2.  Ccanpaiison  of  predate  sphmndal  and  cemstant  radius  tip  geemietries  with  a 
microtip  transmissiem  electron  micrograph  (TEM)  overiay. 

Figure  3.  Prolate  spheroidal  coordinate  syston. 

Figure  4.  Electric  field  vs.  ^  for  a  predate  sidieroidal  tip  with  d=l(X)00  A,  a=10025  A, 
<^=4.7  cV,  and  volt^es  (a)  V=100  volts,  (b)  V=150  vedts,  (c)  V=200  volts. 

Figure  5.  Current  density  vs.  %  for  a  prolate  sidieroidal  tip.  d=10000  A,  a=10025  A, 
<j)f=4,7  cV,  and  voltages  (a)  V=100  volts,  (b)  V=150  volts,  (c)  V=200  volts. 

Figure  6.  Tip  current  vs.  ^  for  a  isolate  sphoxddal  tip.  d=10000  A  a=10025  A, 

4>=4.7  eV,  and  voltages  (a)  V=100  volts,  (b)  V=150  volts,  (c)  V=200  volts. 

Figure  7.  vs.  ^  fw  a  prolate  spheroidal  tip  widi  d=10000  A,  and  selected  values 
of  the  cm  axis  tip  radius  r.  (a)  r=10  A,  (b)  r=50  A,  (c)  r=200  A,  and  (d) 
r=1025  A,  (e)  r=oo. 
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Abstract 

ProMM  variation  inherent  in  the  fiibrication  of  field  entiuer  airaya  leads  to  ntierotip  ndii 
variations  which  can  significanfly  affect  the  device  characteristics  of  vacumn  field 

— ■Hi,  paper  treats  the  tip  laditB  as  a  random  variable  and  exainines  the  effect 
of  tip  Bdins  variation  on  cuirent  emission  wifliin  the  Fowlet-Nordhom  tiieoiy.  Hie 

of  the  sharpest  tips  in  an  stray  is  tpantified  by  nsing  eaunplee  of  an^s 

wMiGamsian  and  Rayleigh  tip  radiiis  probability  densiWfimctioiis  (pdf).  Bfaproposed 

that  tile  K^leigh  pdf  may  be  a  good  model  for  the  real  tip  radhia  variation.  Some 

expected  values  for  tip  ctirient  are  tabulated  fin  a  variety  of  statistical  paranietns. 
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L  Introduction 


Hie  successful  j^cation  of  field  emitter  arrays  (FEA)  as  flat  pand  displays  and  thdr 
potential  application  in  radio  fiequency  (if)  generation  and  amplification  has  spurred 
interest  in  die  yidd-linriting  mechanisms  and  Ae  rdiability  of  tiicse  <kvicc8. 

Reid  emitter  arrays  must  perform  fiieir  function  to  dectiical  specifications,  such  as  die 
total  arr^  current  and  gating  fiequency,  diroughout  their  intoided  lifetimes.  An 
important  issue  for  die  utilization  of  FEAs  is  the  statistical  variability  in  fidmeadon  that 
causes  die  dectiical  parameters  to  vary  fitnn  emitter  to  emitter  [1-4].  These  obser¬ 
vations  show  that  an  of  the  onitters  within  an  array  may  not  function  oc  diat  a  few  of 
die  emitters  are  jHxxJucing  most  of  the  obs^ved  current  mdnle  odiers  are  emitting  at  a 
much  lower  current  level. 

In  a  |MX)peity  fabricated  array  aU  of  the  emitters  are  expected  to  be  fully  functional. 
There  are,  however,  yield  loss  mechanisms  that  limit  tte  total  numba*  of  fully  fac¬ 
tional  devices.  These  loss  mechanisms  typically  fidl  into  flnec  cat^ories.  The  first  two 
are  efaxant  problems  and  random  point  defixts  attributed  to  dust  or  particulate 

fallen  on  the  device  area  during  j&bncatimi.  The  third  yidd  loss  mechanism  is  die  mre 
addressed  in  tins  paper,  namely,  die  parametric  inxxxssing  vaiiatnlity. 
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In  this  paper  we  examine  the  effects  of  statistical  process  variabffity  on  tbe  emission 

current  We  treat  the  tip  radius,  r,  as  a  random  variable  with  a  Normal  (Gaussian)  or 

a  Rayleigh  pdf.  The  pdf  for  the  emission  current  is  subsequently  derived. 

The  results  presented  here  clarify  previous  observations  that  current  emission  comes 
from  only  a  few  tips  in  an  array.  More  impoitandy,  tip  current  expected  values  are 
shown  for  some  selected  arnty  parameters  based  on  the  generally  accepted  Fowler- 
Nordheim  dieory  for  turmeling  current 


n.  Parametric  Process  YmiaWIify 

Ihe  frbrication  of  field  wnission  arrays  (FEA)  is  subject  to  flic  same  process  variation 
inherent  in  all  integrated  circuit  processes  running  under  statistical  process  control. 
Processing  latitude  in  photolithograplty,  wet  chemical  or  plasma  etching,  and  metal  or 
dielectric  deposition  conditions  aU  contribute  to  the  variability  of  geometric  features  on 
the  wafer.  These  process  control  variances  lead  to  FEA  feature  size  variabUity  wthin 
an  array,  from  array  to  array  widiin  a  wafer,  and  from  wafer  to  wafo. 

The  highly  nonlinear  current-voltage  (I-V)  characteristics  of  FEA  devices,  using  eifeer 
die  classic  planar  Fowler-Nordheim  or  the  non-planar  hyperboloidal  tip  models  [5,6], 
magnify  die  effects  of  process  variations  as  seen  in  die  measured  I-V  rdatronship.  One 
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such  exsmplc  is  flic  dependence  of  flie  suifsce  fidd  m  flic  ntfiut  of  curvAtuie.  The 
fifliricdtion  process  steps,  with  their  respective  vsnsnctt  Icsding  to  the  completed 
geomctiy  of  flic  tip  radius,  win  each  contribute  variation  to  flic  final  geometric  sh^ 
of  flie  tips.  Hence,  the  variation  in  tip  radii  for  the  final  stnicture  is  a  function  of  (me 
or  more  independent  variations  due  to  flie  infividual  processing  steps.  The 
independence  of  the  process  variances  can  be  seal  fiom  flic  independence  of  flie  steps 
themselves. 

As  an  exanqile,  consider  that  a  target  feature  size  formed  in  photoresist  depends  (m  flie 
thickness  of  the  photoresist,  flie  exposure  dosage,  and  the  devdopmait  time  med  to 
form  the  desired  image.  Once  flie  image  is  formed  and  flic  ecqiosed  photoresist  is 
removed,  the  image  of  flic  desired  fieatore  is  transfiared  to  the  wafer.  Hw  final 
geometrical  feature,  however,  has  not  been  fcuiued.  Doe  to  variations  in  flie 
photolithographic  steps,  the  im^e  formed  across  the  wafer  surface  exhilrits  variation. 
The  geometrical  feature,  as  an  example  the  gate  opening,  is  fli«i  formed  by  plasma  <x 
wet  chemical  etching.  There  already  exists  a  variatiem  in  flic  {fliotoresist  pattern  doe  to 
photolithographic  processes.  Superin^Kwed  on  top  of  that  photolithogr^faic  variation 
is  the  statistically  independent  vanatitm  due  stflely  to  flie  plasma  or  wet  etching  itsdf 
which  can  originate  from  difiusion  effects  in  the  etchant  chemical,  or  pressure  and 
power  variations  in  the  case  of  plasma  {xocessmg. 
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Tlie  statistical  variatioii  of  a  feature  in  an  array  must  often  be  detennined  by  optical, 
scanning  electron  micrograph  (SEM),  or  transmission  electron  micrograph  (JEM) 
inspections.  These  exhaustive  in-process  inspections  carried  out  during  device 
fabrication  are  time  consuming  and  expensive.  The  most  easily  acquired  statistical^data 
on  arrays  is  obtained  by  using  computer  controlled  testing  of  the  I-V  relationship  for 
completed  devices.  On  dicse  devices  the  overall  effects  of  the  combined  process 
variations  give  rise  to  variations  in  the  measured  I-V  data.  Testing  of  individual  devices 
and  arrays  of  devices  is  perfonned  by  computer  controlled  parametric  testers  which 
record  fee  data  and  can  display  it  in  statistical  form. 

When  seen  from  die  perspective  of  the  conqileted  and  tested  devices,  the  fixmal 
problem  in  tins  case  is  the  relationship  between  the  statistical  variation  in  the  measured 
I-V  data  for  individuany  sanqiled  devices  and  die  variation  in  a  critical  geometric 
feature  within  the  device.  The  formal  analysis  developed  here  can  be  extended  to 
predict  a  statistical  variation  necessary  to  meet  a  total  current  specification  for  a  system. 


HL  Fowler-Nordbeim  Theoiy  for  Mlcrodps 

Rdating  the  dectric  fidd  atdie  surfece  of  a  sharp  lip  to  the  tip  twfius  and  todie  ^Hed 
voltage  has  received  considerable  attention.  Much  of  die  theoretical  work  [7]  has 
adopted  die  idationship  F=pV,  where  F  is  the  dectric  fidd  at  the  emitter  surfece  and 
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V  is  the  appUed  voltage.  The  P  factor,  cMier  comtant  or  (fcpendcnt  <m 
is  used  to  account  for  the  particular  tip  geometry  in  Ae  otherwise  one  dimensional 
Fowler-Nordheim  analysis  [8].  One  theoretical  analysis  [9],  undertaken  to  fit 
experimental  data  based  on  the  Spindt  emitter,  used  a  model  that  invdved  two 
empirically  derived  fectors  to  relate  the  electric  field  to  the  separate  gate  and  collector 
voltages.  Another  effort  in  modeling  [10]  indicated  fliat  the  p  fector  qrproach  may  be 
reasonable  for  tips  with  a  radius  that  is  greater  than  100  angstroms,  but  in  the  case  of 
atomically  sharp  tips  the  P  factor  may  be  an  over-simplification  [8].  A  more  recent 
analysis  [11]  used  a  floating  sphere  model  (FSM)  and  derived  a  flmnula  for  die  p  fv^tor 
based  on  geometrical  considerations.  Dcsirite  the  recent  discussions  regarrfing  flic 
theoretical  necessity  for  die  p  fwtor,  it  is  stfll  widdy  used  as  a  common  basis  for 
con^parison  of  experimental  results. 

In  the  analyses  that  follow,  the  focus  is  on  die  tip-to-tip  radius  variation  in  an  array  and 
how  that  variation  affects  bodi  the  tip  and  aggregate  array  current  For  the  purpose  of 
these  calculations  the  p  factor  will  be  defined  as  p(r)=lflcr.  The  k  term  is  a  geometric 
fector  and  the  r  dependence  provides  a  p  value  diat  increases  with  decreasing  ra<fius  and 
yidds  greater  fidd  enhancement  than  a  constant  P  value.  The  k  tom  will  be  assumed 
to  have  a  constant  value.  The  dectric  fidd  and  radius  rdadonship  can  now  be  used  in 
the  FowlCT-Nordheim  equation  to  calculate  the  current  doisity  for  a  dp.  In  diose  cases 
where  p  is  defined  differently  from  p(r)=iykr,  the  cxpUcit  dependence  of  p  on  the 
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geometrical  feature  size,  r,  is  incoiporated  in  the  expression  for  foe  wuueat  density  j(r). 
The  derivation  of  foe  cunent  probabiUty  density  function  then  proceeds  in  foe  same 

manner. 


Tbc  most  general  form  of  foe  Fowler-Nordhcim  equation  is  [12]: 


a) 


In  fois  expression  F  is  foe  dectric  fidd  (V/anX  ♦  »  foe  material  woric  function  (eV), 
h  is  Plank's  constant  (ag-secX  e  is  foe  dectron  chaige  (stat-coulX  and  m  is  foe  dectron 

mass  (gm). 

The  V  function  is  dependent  on  an  involved  conqwsition  of  elliptic  int^rals  and  varies 
from  1  to  0  as  dectric  fidd  strengfo  increases.  The  t*  function  is  a  slowly  varying 
function  of  foe  applied  fidd  and  is  approximated  as  1.1  [12].  In  ordCT  to  avoid 
calculations  wifo  foe  dfiptic  integrab,  foe  function  v  is  approximaled  by  [13]: 

=  QJ95-jiFf  j{F)=3.T9xlO'*^  (2) 
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Using  this  approximation  and  substituting  the  values  for  t*  and  die  constants,  j(F)  can 
be  expressed  in  an^cm*  as: 


m' 


IMxlO^F^ 


03) 


This  form  of  die  equation  uses  electric  fidd  values  in  volt/cm.  Replacing  die  electric 
field  in  fiivor  of  the  applied  voltage  V  and  transforming  to  die  radfius  of  curvature  in 
angstroms  as  die  independent  variable  r,  die  current  density  expression  becomes. 


where  the  constants  A  and  B  for  a  given  V  are  defined  asi 


U4ciO^ _ 9.811 

fPCP 


B=^ 
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0.649^^ 

V 


IV.  Caiciihition  of  Total  Current 

The  Fowler-Nordheim  equation  developed  above  provides  for  tbe  current  density  j(r)  in 
ampiJar?  for  each  individual  tip  widi  its  particular  radiiis.  The  cuiient  fiooi  a  microtip 
can  be  found  by  multiplyiiig  the  current  dmsity  j(r)  by  die  eSeedve  emitting  area  for 
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diatniiiis.  Theto«dcunaitofa»I^it»»»nni«lionofl^ciin«nlfcr«nofaio 

individual  Spa.  Alftough  the  FowlavNordheim  and  noniianar  emiaaion  Iheoriea 
expreaa  die  emiaaion  cunent  aa  a  crnieni  denaily,  the  total  emiaaion  cutrent  from  an 

individual  Spot  an  anayia  the  parameter  of  inteiest  Aa  auch,  we  rewrite  the  cunent 

tosily  aa  a  Sp  current  using  an  appropriately  chosen  expreaaion  fcr  die  emiaaion  area. 

There  ia  no  general  agreement  on  the  calculaSon  of  effecSve  enritting  area  for  a 
nucroSp.  A  recent  theoreScal  atudy  [9]  uaed  the  modd  of  a  apherical  Sp  on  die  end  of 
a  cone  and  calculated  the  area  of  emission  widi  a  dependence  on  die  emitter  Sp  half 

angle. 

The  approach  shown  in  the  calcrdaSons  diat  Mow  is  to  adopt  an  emission  area  for  the 
Spa  that  can  be  viewed  aa  Ireing  dependent  on  the  micrutMcaSon  procedure  used  to 

create  die  Sp.  hi  parScular,  our  experimentaSon  on  atomically  sharp  silicoo  emitten 
made  use  of  a  technique  to  isotropically  etch  silicon  while  using  a  sOicon  cBoxide  Tiat* 

asamask.  This  technique  creates  concave  sides  for  the  nScroSp  structure,  and  these 

sides  would  ideally  he  perpendicular  tothefineSpareaatlhentoat  appa  part  of  the 

sidewaa.  Thus,  the  enSssion  area  drat  is  used  in  the  ftBowingcalcolaSoiisia  the  half- 

spherical  ahdl  aurfice  deswibed  by  a(r)=2ia>,  where  r  is  expressed  in  cenSmeteis. 

The  Sp  current  in  amperes  can  then  be  expressed  as  ifrHfrMr). 
following  expressim: 
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^)=2xlO~^%A  «xp(-flr) 


(0 


The  radius  r  is  in  angstroms.  In  most  geometrical  cases  die  effective  emitting  area 
should  scale  in  proportion  to  die  square  of  the  radius,  certainly  to  a  good  approximation 
if  not  exactly.  'Iheiefore,  diis  sinqiliricadon  in  die  final  expression  for  i(r)  should  exist 
in  most  tip  geometries. 


V.  Cment  Probability  Density  Function  Derivation 

The  tip  raifius  can  be  viewed  as  a  conlinuoos  random  variable  in  a  given  range.  This 
random  variable  is  described  by  a  probability  density  function  (pdf)  over  diat  range. 
If  3JJ  <>ppropri(*te  tip  radius  pdf  can  be  determined,  dien  sto^iastic  methods  [14]  can  be 
used  to  analyze  field  emission  parameters  in  light  of  the  Fowler^Nordheim  equation  and 
non-planar  emission  theories  [5,6]. 

We  considered  two  airtys  of  microtips  in  which  the  ra^fius  is  modded  as  a  continuons 
random  variable^  die  first  distributed  Gaussian  second  distributed  Raleigh 

hR(r)  as  shown  in  Table  1.  In  all  cases  the  rai^e  of  radii  values  is  100-300  A 
(angstroms). 
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To  find  die  cuncnl  probabifity  doisiv  fimOioiB  m  pdf  »nd  »» 

dK  Rayidgh  pdf)  deacribiiig  cmrent  for  «  single  lip,  flie  roof  (r,)  for  the  twfius  r  and  die 

derivative  of  the  current  must  be  found  [15]  from  equation  (6): 


AiQ-')=-2xlO''*%AB  exfi-Br)  (8) 

dr 


TTie  corresponding  current  probability  density  functions  are  shown  in  Table  1.  These 
functions  are  plotted  in  Figure  1.  The  total  probabiUty  area  under  each  of  the  curves 
has  been  calculated  and  shown  to  be  unity.  Since  the  abscissa  current  values  are  in 

amperes  and  cover  a  range  of  only  about  0.001  anipere,  Ae  plotted  ordinate  values  can 

be  exceedingly  large,  hi  the  continuous  random  variable  case  die  point  probal^ties  at 
a  particular  entrent  value  have  no  meaning;  only  in  a  range  of  current  values  can  a 
determination  of  probabiUty  be  made.  The  expected  values  for  tip  current  were 
computed  and  are  also  shown  in  Rgure  1. 
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Using  a  »«TnnM-  procedure  for  finding  fimcticaa  of  •  nodom  variaUe,  die  curratf 
probability  density  functions  describing  the  total  current  fiom  an  array  are  derived.  The 
resultant  total  current  probability  density  functions  minor  the  tip  current  probability 
density  functions  with  the  substitution  of  die  total  array  current  variable  for  the  tip 
current  variable.  The  expected  values  for  the  total  array  current  are  those  for  the  tip 
current  multiplied  by  die  number  of  tips  in  die  array. 

Choosir^  between  probabihty  density  functions  fiir  modeling  the  array  current 
performance  is  inqiortant  Since  the  actual  distribution  is  ififficult  and  expensive  to 
obtain,  it  often  ranains  unknowrL  Howeva:,  some  insight  firan  expaimentation  may 
be  used  to  identify  a  likely  pdf  for  a  givoi  fidiiicatian  process. 

We  propose  the  use  of  die  Rayleigh  distribution  to  lepresoit  die  tip  radius  variation. 
Consider  that  etch  rates  also  vary  as  a  function  of  crystal  {darK  or  lattice  dtiection.  For 
example,  in  a  simple  cubic  lattice  die  directions  x,  y,  and  z  could  experience  different 
etch  rates.  If  diese  particular  coordinates  at  the  tip  surfecc  are  considered  Gaussian 
distributed  random  variables  due  to  etching  variations  along  dior  respective  directions, 
then  die  random  variable  representing  die  radius  would  be: 
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With  tins  dcMticm.  the  nmdom  variable  R  would  be  Rjyld^  'Hw 

same  ^jproach  can  be  generalized  to  include  any  number  of  lattice  directions  (he^  for 
the  diamond  lattice)  and  the  resulting  tip  radius  pdf  would  still  be  shown  to  be  Raylagh. 


The  Rayleigh  dstribution  is  shown  plotted  against  Ae  Gaussian  distribution  in  Figure 
2,  The  Raleigh  parametes  have  been  set  to  align  closdy  witti  the  Gaussian 
distribution  so  diat  a  con^arison  of  current  calculations  can  be  made.  The  conjecture 
presented  here  is  tfiat  while  die  Gaussian  tip  radius  pdf  is  of  interest,  die  Rayleigh  pdf 
is  probably  a  more  realistic  choice  in  moddling  the  field  emission  current  from  an  array. 
In  a  f?Kr;^afirtn  process,  die  photoUthography  and  etching  steps  are  designed  to  achieve 
a  targeted  tip  rarfius  value.  As  the  process  wet  or  plasma  etching)  is  carried  out 
to  achieve  the  limits  of  its  capabilily  0.c,  sharpness^  die  process  variation  will  be 
exhibited  in  die  laoccss  parameter  diat  represents  d^radadon  fitmi  its  liimts  ^e., 
decreased  tip  sharpness).  The  resultant  tip  radius  pdf  would  exhibit  a  cutoff  at  the  small 
radii  end  of  die  pdf  radier  than  the  tail  that  characterizes  the  Gaussian  distribution. 


VIL  Summaiy  of  RcstiKs 

The  variation  in  tip  radius  inherent  in  the  fiJnication  of  a  field  emission  array  is  a  result 
of  the  statistical  nature  of  die  &brication  process.  The  resultant  I-V  characteristic  is 
strongly  dependent  on  the  sharpest  tips  in  the  array.  To  ^lantify  die  performance  of  die 
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airay,  Ac  tip  radius  is  viewed  as  a  random  variable  and  stochastic  concepts  are  applied. 
In  Ac  case  of  a  Gaussian  distributed  tip  radius,  Ac  array  current  has  a  lognormal 

distribution. 

It  can  be  observed  that  Ac  tips  wiA  a  small  radius  have  a  dominating  effect  relative  to 
Ac  larger  radius  tips.  For  Ac  Gaussian  array,  Ac  tips  wiA  radii  between  100-150  A 
contribute  about  1%  of  Ae  current  because  Acre  are  so  few  tips  in  that  range  for  Ac 
chosen  standard  deviation  of  15  A;  Ae  tips  150-200  A  contribute  85%,  Ac  tips  200-250 
A  contribute  14%,  and  Ac  tips  250-300  A  contribute  a  negligible  current  percentage. 

It  has  been  proposed  that,  at  least  in  some  cases,  Ae  variation  of  Ae  tip  ra<&B  may  be 
more  closely  modeled  by  a  Rayldgh  distribution  raAer  than  a  Gaussian  (fistribution. 
The  Rayleigh  tip  radius  pdf  exanqile  results  in  an  array  current  Aat  is  about  60  percent 
of  Ac  current  calculated  from  a  similar  array  wiA  a  Gaussian  Astributed  tip  taAus. 
Additional  calculations  of  expected  tip  current  for  various  parameters  are  listed  in  Table 

2.  Hic  particular  Rayleigh  pdf  parameters  listed  to  the  right  of  Ac  Gaussian  parameters 

have  been  determined  to  match  Ae  two  Astribution  functions  as  shown  in  Figure  2. 
Comparison  of  Ae  horizontal  pairs  of  expected  tip  currents  inAcates  Aat  Ae  Rayleigh 
expected  current  falls  below  Ac  Gaussian  expected  currait,  and  to  a  greater  degree  as 
Ae  Astribution  spread  or  standard  deviation  is  increased. 
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Hie  tabulated  values  of  expected  current  are  plotted  in  Hgure  3  to  hig^iligbt  the  ratio 
of  the  standard  deviation  to  Ac  mean.  The  Rayleigh  curves  are  consistently  below  Aeir 
respective,  paired  Gaussian  curves.  The  upward  trending  of  tiie  Gaussian  curve  (a)  in 
Figure  3  is  reflective  of  Ae  feet  tiiat  Ac  pdf  functions  wiA  Ae  ratio  o/ji=0.1  have  Acir 
Gaussian  tails  substantially  wiAin  Ae  limits  of  tip  radii  used  to  add  iq>  current 
contributions.  The  increased  standard  deviations  of  Ae  oAer  Gaussian  curves 
designated  as  (b)  and  (c)  (i.e,  o/n=0.2,  o/n=03)  resulted  m  a  substantial  portion  of  Ae 
pdf  Gaussian  felling  outside  Ac  fimits  of  tip  radii  used  to  add  current 
cerntributions. 

The  practical  i"atter  of  examining  Ae  ndcrotqis  of  an  array  by  sanqiling  to  determined 
wheAer  Ae  t^  rarfius  is  distributed  eitiiCT  Gaussian  or  Rayleigh  is  indeed  difticult  As 
shown  in  Figure  2,  such  a  task  is  a  detenninatirm  of  wheAer  Ae  Gaussian  tail  is  present 
or  not  An  alternative  for  making  such  a  pdf  determination  is  to  sample  Ac  microtips 
for  Ae  purpose  of  determining  a  mean  value  radius  and  a  standard  deviatioiL  This 
statistical  information  can  tiien  be  used  to  modd  Ae  microtip  array  as  demonstrated 
hodn,  and  tiic  modeled  current  can  be  conqiared  wiA  actual  current  in  maku^  a  final 
determination  of  tiie  tip  radius  pdf. 

It  is  important  to  note  that  Ae  above  analysis  si^ests  tiiat  it  may  be  difScult  to 
accurately  predict  Ae  performance  of  an  array-^^pe  emission  device,  based  on  tiie 
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measured  enrissicm  characteristics  of  a  sin^  We  have  shown  that  the  total  current 

of  such  a  device  depends  in  a  nonlinear  fishion  <mi  the  statistical  variation  of  the  tip 
radii  which  will  inevitably  occur  in  a  &brication  process,  and  it  is  acc«iluated  by  using 
nricron-sized  procedures  to  form  nanometer-sized  structures.  However,  the  understanding 
and  quantification  of  the  statistical  effects  associated  with  tip  radius  variation  can  assist 
in  developing  a  device  testing  strategy  to  support  a  j&brication  specification. 
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CAPTION  PAGE 


Figure  1.  Tip  current  probability  densi^  functions  for  Gaussian  and  Rayleigh 
distributions  of  tip  radii  between  100-300  A,  with  V=150  volts,  k=1.6,  and  eV. 
(a-Gaussian,  b-Rayleigh).  Calculated  expected  values:  Gaussian  2,16  pA  and  Rayleigh 

136  pA. 

Fgure  2.  Conqiariscm  of  a  Gaussian  pdf  to  a  Ra^dgh  pdf  for  the  tip  range  10^ 
300  A.  The  Gaussian  parameters  are  mean  p=200  A  and  standard  deviation  o=15  A. 
The  Rayleigh  parameters  are  y=175  A  and  X=23  A.  (soUd-Gaussian,  dashed-Raylcigh). 

Fgure  3.  Expected  tip  current  vs.  tip  radius  for  various  o/p  ratios.  (a-Gaussian  with 
o/p=0.1,  b-Ganssian  with  o/p=03,  c-Gaussian  with  o/p=03,  d-Rayld^  similar  to  curve 
a,  e-Rayleigh  Mmilar  to  curve  b,  f-Raylcigh  similar  to  curve  c). 

Table  1.  Conqiarison  of  tip  current  probability  density  functions  based  on  Gaussian  and 
Rayleigh  tip  radius  probability  density  functions.  O  is  the  Heaviside  unit  step  function 
in  die  Rayleigh  pdf. 

Table  2.  Tip  Current  Expected  Value  Exanqiles. 
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Appendix  D 
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I.  INTRODUCTION 

Experiments  on  frequency  mixing  [1*3]  and  tunneling  time  measurements  [4-6] 
have  stirred  renewed  interest  in  the  problem  of  tunneling  through  time  dependent  barriers 
in  a  vacuum  diode  configuration  such  as  the  STM  junction.  When  the  laser  field  interacts 
with  the  STM.  the  tip  and  base  act  as  an  antenna  system  for  detection,  harmonic 
generation  and  mixing  [7].  The  resulting  nonlinear  l-V  characteristics  of  the  STM  diode 
produce  both  a  rectified  DC  and  AC  tunneling  voltage  bias  [7].  Earlier,  the  metal-oxide- 
metal  point-contact  diode  with  the  same  configuration  was  used  to  determine  absolute 
frequencies  up  to  the  visible  (-10'®  Hz)  and  as  a  consequence  of  these  precise 
measurements,  the  speed  of  light  is  redefined  as  an  exact  SI  constant  of  299792458  m/s 
[8].  Lucas  et  al.  [9]  and  Krieger  and  co-workers  [1-3]  have  claimed  that  the  STM 
configuration  is  superior  to  the  point-contact  diode  because  of  its  controllable  gap  width 
and  electron  transport  in  vacuum.  This  latter  is  easier  to  treat  theoretically  than  transport 
through  an  oxide  layer.  However,  at  present  the  frequency  response  of  the  STM  junction 
has  only  been  measured  to  a  few  tens  of  gigahertz  (-10'®  Hz).  There  are  few  theoretical 
treatments  of  electron  transport  through  oscillating  barriers  [10-14].  Here  we  describe 
briefly  the  formulation  of  a  theory  and  present  calculations  of  the  inelastic  photo-assisted 
tunneling  through  a  planar  junction  with  both  a  static  DC  voltage  and  a  AC  voltage  of 
angular  frequency  w.  One  purpose  of  this  paper  is  to  ultimately  explain  the  experimental 
data  of  Nguyen  et  al.[4]  for  the  determination  of  the  tunneling  time  in  a  laser  irradiated 
STM  junction. 

In  1987  Cutler  et  al.  [15]  proposed  an  experiment  to  estimate  the  time  of  transit  of 
an  electron  tunnefing  through  a  quantum  mechanical  barrier  by  measuring  the  DC  current 
as  a  function  of  gap  spacing  (for  a  fixed  laser  frequency)  or  as  a  function  of  frequency 
(for  a  fixed  gap  spacing)  in  an  irradiated  STM  junction.  If  a  linearly  polarized  laser  beam 
is  focused  on  the  STM  junction  with  the  shank  of  the  tip  assumed  to  be  many 
wavelengths  longer  than  the  radiation,  then  the  tip  acts  as  an  antenna  to  couple  the 
radiation  and  the  junction.  The  current  at  the  laser  frequency  will  be  propagated  with 
sntall  attenuation  over  an  appredable  length  and  establish  a  time  dependent  osdllatirig 
potential  between  the  tip  and  base.  Because  of  the  geometiic  asymmetry  of  the  junction, 
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a  DC  component  is  generated  of  the  form 


2"(2n)!dV^' 


(1) 


where  the  second  term  arises  from  the  non-linearity  of  the  junction  [4].  To  estimate  the 
tunneling  time  they  assume  a  simple  pseudo-kinetic  description  of  ballistic  transport,  for 
the  case  of  fixed  gap  distance  at  a  constant  frequency.  If  the  junction  is  irradiated  by  a 
laser  of  a  fixed  frequency  where  the  polarization  is  parallel  to  the  axis  of  the  tip,  the 
electron  transfer  process  is  enhanced  during  half  of  the  period  when  the  electric  field 
vector  of  the  laser  beam  Is  accelerating.  During  the  second  half  of  the  period  when  the 
laser  field  is  reversed,  the  barrier  will  Increase  and  transport  will  be  diminished  [4,9].  If 
the  gap  width  is  greater  than  some  critical  value  s^,  the  electron  will  not  have  "enough 
time"  to  traverse  the  barrier  and  be  detected  before  the  field  is  reversed.  The  tunneling 
time  can  be  deduced  from  [4.16] 


V 


(2) 


where  It  is  found  that  vssvp,  the  Fermi  velocity  of  the  electrons  in  the  metallic  electrodes 
[4,16].  The  argument  for  estimating  the  tunneling  time  is  similar  for  the  case  of  a  fixed 
gap  spacing  with  a  variable  frequency.  If  the  frequency  changes  too  rapidly  then  the 
electron  will  not  have  "enough  time"  to  traverse  the  barrier  before  the  field  reverses 
direction. 

In  1989,  Nguyen  et  al.  [4]  used  the  dynamical  response  of  an  STM  to  measure  an 
operational  tunneling  time.  In  this  method,  a  natural  time  scale  is  provided  by  the 
frequency  of  the  laser  that  constitutes  an  integral  part  of  the  measuring  apparatus.  The 
laser  causes  the  tunneling  and  at  the  same  time  determines  the  "speed"  of  the  tunneling 
process.  In  this  experiment,  a  linearly  polarized  Nd-YAG  laser  of  wavelength  i=1 .06Mm 
was  focused  on  an  STM  Junction  consisting  of  a  W-tip  (of  radius  of  curvature  -2nm  and 
a  Si  base  in  UHV  (-1(r’®lorr).  The  experiment  was  performed  with  no  applied  DC  bias. 
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The  rectified  current  was  measured  from  tip  to  sample.  From  their  results,  using  Eq.  (2), 
they  deduced  a  tunneling  time  of  about  1 .8  fs.  This  is  in  agreement  with  the  theoretical 
estimates  of  Hartman  [17^,  Buttiker  et  al.  [11]  and  Huang  et  al.  [16].  A  somewhat 
different  version  of  this  experiment  by  Moller  et  al.  [18]  yielded  similar  results  but  the 
interpretation  of  the  tunneling  time  is  still  in  question.  Other  attempts  have  been  made  to 
measure  the  "traversal  time"  [5,6].  Most  related  Is  the  experiment  of  Gueret  et  al.[5]  who 
indirectly  measured  the  traversal  time  in  a  heterostructure  tunneling  barrier  using  the 
orbital  analog  of  a  Larmor  spin  precessional  dock.  Their  value  of  t  for  a  barrier  of  40-50 
nm  in  width  scales  with  the  results  of  Nguyen  et  al.  for  a  barrier  of  1-2  nm  [4]. 

In  this  paper  we  first  describe  briefly  the  formulation  of  a  theory  of  the  inelastic 
mediated  photo-assisted  tunneling  for  a  planar  bimetallic  junction  with  both  a  static  DC 
voltage  and  a  time  dependent  bias  of  angular  frequency  o>.  A  more  detailed  derivation 
and  analysis  will  be  published  elsewhere  [14].  We  then  present  results  for  a  trapezoidal 
bam'er  surmounted  by  an  oscillating  potential  of  the  same  frequency  range  as  used  in  the 
tunneling  time  experiments  [4]. 

In  Section  II  an  expression  for  the  current  Is  derived  using  a  kinetic  formulation  of 
the  transfer  hamiltonian  method.  Although  the  current  is  initially  calculated  within  the  free 
electron  model,  the  kinetic  formulation  allows  for  the  indusion  of  band  structure  effects 
through  the  density  of  states  factors.  In  section  III  the  formalism  is  applied  to  a  model 
planar  MVM  diode  with  a  trapezoidal  barrier  (due  to  an  external  static  field)  surmounted 
by  a  time  dependent  oscillating  potential.  Although  the  wave  function  can  be  calculated 
exactly  to  infinite  order,  we  have  initially  calculated  the  current  to  first  order  in  order  to 
focus  on  the  role  of  single  photon  mediated  inelastic  tunneling.  The  transition  probability 
is  calculated  using  time  dependent  perturbation  theory.  The  photo-assisted  inelastic 
tunneling  current  was  obtained  as  a  function  of  frequency,  fields,  barrier  parameters  and 
material  asymmetry.  The  results  presented  in  Section  IV  indicate  that  for  the  tunneling 
biases  used  in  STM  experiments,  the  contribution  of  inelastic  mediated  photo-assisted 
tunneling  rises  rapidly  for  frequencies  a  10’*  Hz  (I-©-,  visible).  The  results  are  applied  to 
an  analysis  of  the  tunneling  time  experiments  of  Nguyen  et  aL  [4].  The  conclusions 
suggest  some  constraints  but  do  not  invaTidate  the  experimental  procedure  used  for 
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extracting  tunneling  times  in  an  irradiated  tunneling  junction.  This  is  discussed  in  Section 
V. 


II.  TUNNELING  CURRENT  FORMALISM 
The  tunneling  current  is  calculated  using  a  kinetic  formulation  of  the  transfer 
hamittonian  method.  Time  dependent  perturbation  theory  is  used  to  determine  the 
perturbed  barrier  wave  functions.  The  harmonic  perturbation  is  of  the  form, 

V(t)*eV^coso)t,  (3) 

where  Is  the  amplitude  of  the  applied  oscillating  field  and  o  is  the  angular  frequency. 

The  inelastic  (and  elastic)  tunneling  enhancement  due  to  this  perturbation  is  viewed 
as  a  transition  from  an  electron  energy  state  E^,  (electrode  to  the  left  of  the  barrier)  to  an 
energy  state  E^,=E^±m>i&>  (electrode  to  the  right  of  the  barrier),  where  m  is  an  integer. 
For  weak  coupling  of  the  two  electrodes,  the  wave  function  inside  the  barrier  region  can 
be  written  as. 


T(r.t)=a^<t>^(r)e“+j;  b  Jt)(D^(t)e~ 

m 


(4) 


where  a(t)  and  b(t)  are  time  dependent  coefficients.  The  functions  0^(r)  and  <J>^(r)  are 
the  spatial  solutions  of  the  unperturbed  left  and  right  hamiltonians  H,  and  respectively. 
These  barrier  wave  functions  have  decaying  exponential  like  tails  in  the  bamer  region 
whose  overlap  is  responsible  for  the  tunneling. 

Using  time  dependent  perturbation  theory  and  the  initial  conditions  a(0)=1  and 
b{0)=0,  we  obtain  to  first  order  [19], 


(5) 


Here  transition  rate  from  the  initial  state  (0  to  the  final  state  (f)  is 
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We  assume  that  the  final  states  are  dense  enough  so  that  the  total  transition  probability 
in  Eq.(6)  can  be  expressed  in  terms  of  an  integration.  If  the  parallel  component  of  the 
energy  is  conserved  during  tunneling,  the  transition  rate  from  state  to  E,  can  be  written 

R,>IjdE„p(E„)|b(t)|»(1-f(EJ),  (7) 

where  the  x-direction  is  the  normal  direction  and  the  Fermi  factor  has  been  introduced 
to  account  for  occupancy  of  the  final  state.  Since  the  transmitted  probability  current  Is 
equal  to  the  tunneling  rate  the  transmission  coefficient  or  tunneling  probability  is 


(8) 


where  Rq  is  the  probability  current  inddent  on  the  barrier.  For  a  free  electron  model  and 
box  normalization, 

Ro=^|a„l'.  (9) 

m 

where  the  coeffident  of  the  inddent  plane  wave  is  £^,. 

The  tunneling  current  within  the  kinetic  formulation  is  the  product  of  the  arrival  rate, 
density  of  Initial  states,  occupation  factor  for  initial  state  and  tunneling  probability 
integrated  over  all  initial  states,  E,, 


J»=Jp(Ej)v(E^f(Ei)Dj^Ej,  (10) 

where  v  is  the  arrival  rate  of  electrons  on  the  barrier  and  p{E,,EJ  is  the  density  of  states 
for  Inddent  electrons  having  energy  between  E,  and  E,+dE,  and  normal  energy  between 
E^  and  E^+dE^. 
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Within  the  free  electron  approximation  the  transmission  coefficient  Is  [20] 


vJa  .n 


sin*o!t/2  sin*<i)^t/2 


2 


v,|a„|- 

where  and 

|MJ*s|<(DJeV^j<t>,>|2 

for  a  harmonic  perturbation.  In  the  limit  of  long  time, 


2 


+cross  terms 


(11) 


(12) 


1  sin2(E^V2)^) 

lim - = - 

t—t  (E.t/2h)2 


=2rrt^{EJ 


(13) 


The  delta  function  limit  yields  the  result  that  the  most  probable  transitions  involve  the 
exchange  of  a  single  photon.  The  short  time  limit  Is  proportional  to  t  indicating  that  the 
tunneling  current  initially  increases  with  time  [19]. 


III.  TUNNEUNG  CURRENT  FOR  A  TRAPEZOIDAL  BARRIER  IN  FE  APPROXIMATION 
In  this  section  we  calculate  the  current  for  the  case  of  a  trapezoidal  barrier  within 
the  free  electron  approximation.  The  hamiltonian  for  the  left  electrode  is 

2 

H,=-— v2+[Vo-Fx]0(x-s),  (14) 

'2m 

and  the  left  electrode  hamiltonian  is 

2 

H,=-^v“*tv,-Fx)e(s-x)*[C,-C,-eVje(x-s).  05) 

Here  G(x)  is  the  step  function,  F  is  the  applied  static  electric  field,  are  the  Fermi 
energies  of  the  left  and  right  electrodes,  and  ’s’  Is  the  barrier  gap  width.  The  zero  of 
energy  is  taken  to  be  tiie  bottom  of  the  left  hand  electrode.  The  solutions  for  the 
unperturbed  hamiltonians  are 
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where  A  is  the  normalized  area  of  the  planar  contact  and  k,  is  the  transverse  component 
of  the  wave  vector.  The  solutions  for  are  given  by  a  linear  combination  of  plane 
waves  in  the  electrodes  and  by  Airy  functions  of  the  first  and  second  kind  in  the  barrier 
region  [20,21].  Using  boundary  matching  at  the  two  interfaces,  x  can  be  determined 
exactly  in  each  region  [20,21]. 

The  total  tunneling  current  is  [20], 


f  ^ciE„|dE^JdE,f,(E^[1  -f(E,-E^*E„*eVJ] 
xp{E,,EJp(E„»C,-{,*eVJ[i(E„-E^->,u)l|MJ» 


00  oe 


1 

|b3,l 


-|dE^  I  ^dE„^dE,f/E,.CrC.-eVJ[1-f(E,-E„.EJ 


(17) 


2^ 

') 


The  subscript  1  refers  to  the  left  hand  electrode,  2  to  the  barrier  region,  and  3  refers  to 
the  right  hand  electrode. 


IV.  RESULTS 

To  determine  the  frequency  response  of  a  model  W-Au  junction  we  calculated  the 
tunneling  current  using  Eq.(17).  In  Rg.  1  we  exhibit  the  results  of  the  tunneling  current 
as  a  function  of  angular  frequency  with  a  static  field  F^=1.5x10®V/m  and  an  oscillating 
field  of  magnitude  F^=5x10V/m.  Calculations  were  done  for  two  different  gap  widths, 
s=0.8nm  and  s=1.0nm.  For  <i)<5x10’Vad/s,  the  current  dentsity  is  independent  of 
frequency;  for  o>10'^rad/s,  there  is  a  rapid  increase  in  the  current  density,  i.e.,  about  a 
tenfold  increase  in  J  as  o  goes  from  10’®  to  3x10’®rad/s.  This  behavior  Is  exhibited  fiar 
both  gap  widths  and  other  model  calculations  by  the  authors  give  the  same  systematics 
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for  different  values  of  the  static  field.  The  calculated  results  are  in  good  qualitative 
agreement  with  the  analysis  of  Buttiker  and  Landauer  [11]  who  calculated  the 
transmission  coefficient  for  time  dependent  tunneling  in  an  oscillating  barrier 
superimposed  on  a  one-dimensional  rectangular  barrier.  Their  calculated  transmission 
coefficient  for  the  first  two  side  bands  also  indicates  that  at  low  frequency  the  probability 
of  tunneling  is  independent  of  frequency.  However,  at  higher  frequencies,  the 
transmission  probability  increases  rapidly.  This  is  consistent  with  the  results  in  Fig.  1. 
This  occurs  because  a  particle  which  absorbs  a  quantum  and  has  energy 
tunnels  through  the  barrier  more  easily  than  particles  with  energy  E  or  E-l>o.  For  angular 
frequencies  sufficiently  greater  than  -5x1 0^^  rad/s,  electrons  can  gain  enough  energy  so 
that  pure  photoemission  occurs. 

We  have  also  examined  the  dependence  of  the  tunneling  current  density  on  the 
gap  spadng  with  and  as  parameters.  In  Rg.  2  we  exhibit  the  results  for 
<D=1.8x10^*rad/s  (approximately  the  value  for  the  Nd  YAG  laser)  and  for  a  low  frequency 
o>=5x10^Vad/s.  Both  the  static  and  oscillating  field  are  lOmV,  values  typical  to  some  STM 
experiments.  For  a  constant  V„  the  tunneling  current  density  falls  off  approximately 
exponentially  as  s  Increases,  which  mimics  elastic  tunneling. 

Using  Eq.(17)  we  calculated  J  as  a  function  of  ’s’  at  two  different  frequencies  but 
with  a  constant  F^=1.5x10®V/m  and  two  values  of  F^=1.5x10®V/m  and  F^=5x10^  The 
results  are  plotted  in  Flg.3.  The  curves  exhibit  similar  behavior  with  a  maximum  and  rapid 
fall-off  vrith  gap  distance.  This  results  from  ohmic  transport  for  small  gap  spacings  and 
then  conventional  tunneling  as  the  gap  distance  increases.  The  curves  exhibited  In  this 
figure  are  almost  identical  quantativeiy  to  the  low  frequency  results  of  Levy  Yeyati  and 
Rores  [13].  Calculations  show  that  the  peak  position  and  fall-off  is  relatively  insensitive 
to  frequency  even  up  to  frequencies  near  the  visible.  The  curves  exhibit  a  cut-off  in  the 
rectified  current  as  a  function  of  gap  width  for  a  fixed  a  and  F^  as  was  predicted  by 
Cutler  6t  al.  [15].  However,  the  cutoff  does  not  have  the  expected  dependence  on 
frequency.  One  of  the  reasons  may  be  that  only  a  field  can  reverse  the  direction 
of  the  electron’s  motion  and  lead  to  the  expected  decrease  in  Jqq.  In  these  first 
calculations  the  static  field  is  always  equal  to  or  greater  than  F^.  In  addition,  even  the 
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small  work  function  differences  of  the  W-Au  junction  (-0.2-0.3eV)  correspond  to  effective 
fields  of  about  2-3x1 0®V/m  for  a  gap  width  of  1  nm.  Hence,  the  static  field  always 
dominates. 


CONCLUSIONS 

In  the  present  work  we  have  formulated  a  theory  of  inelastically  mediated  photo- 
assisted  emission  for  a  planar  bimetallic  junction  with  both  a  static  dc  voltage  and  a  time 
dependent  bias.  An  expression  for  the  current  is  derived  using  kinetic  theory  and  the 
transfer  hamiltonian  method.  The  wave  functions  are  determined  from  time  dependent 
perturbation  theory  and  are  calculated  recursively  to  infinite  order.  The  photo-assisted 
inelastic  tunneling  current  was  calculated  as  a  function  of  frequency,  fields,  barrier 
parameters  and  material  asymmetry. 

Some  important  conclusions  emerge  from  these  calculations  which  impose 
limitations  on  but  do  not  invalidate  the  experimental  procedure  to  determine  a  tunneling 
time  using  a  laser  irradiated  STM  junction.  To  make  the  experiment  feasible  it  is 
necessary  that  the  rectified  DC  current  must  reach  a  maximum  and  then  decrease  as  a 
function  of  frequency  for  fixed  gap  spacing  or  as  a  function  of  gap  spacing  for  a  fixed 
frequency.  In  the  first  method,  having  chosen  some  appropriate  gap  spacing,  the  critical 
frequency  at  which  Ijjc  decreases  to  some  negligible  value  must  be  such  that  the 
electron’s  "motion"  can  be  reversed  by  the  oscillating  field  before  it  is  detected  at  the 
anode.  This  requires  that  the  effective  static  field,  consisting  of  the  applied  static  bias 
and/or  the  field  due  to  the  contact  potential  difference  is  less  than  the  amplitude  of  the 
laser  induced  oscillating  field.  In  the  alternative  procedure,  the  limitations  are  less 
stringent  and  only  the  condition  that  (VJ  <F^(V  J  needs  to  be  satisfied.  A  careful  and 
systematic  study  will  have  to  be  done  before  definitive  criteria  can  be  established  to 
define  the  exact  limitations  on  the  experimental  procedure. 
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LIST  OF  FIGURES 

Figure  1  The  tunneling  current  density  as  a  function  of  angular  frequency  for  gap 
widths  s=1nm  (solid  line)  and  s=0.8nm  (dotted  line).  The  average  static  field  is 
1.5x10®V/m  and  the  amplitude  of  the  oscillating  field  is  5x10*V/m. 

Rgure  2  The  tunneling  current  density  as  a  function  of  gap  width  with  angular 
frequency  as  a  parameter.  The  dotted  line  indicates  the  values  for  Q=1.8x10’®rad/s  and 
the  solid  line  is  for  o)=5x10’^rad/s.  The  external  bias  voltage  is  10  mV  and  the  amplitude 
of  the  oscillating  field  Is  10®V/m. 

Figure  3  The  tunneling  current  density  as  a  function  of  gap  width  with  amplitude  and 
angular  frequency  as  parameters.  The  average  static  field  is  1 .5x1 0®V/m.  The  applied 
fields  F^=1.5x10®V/m  ( — )  and  F„=5x10^V/m  for  «=1.8x10’®rad/s  ( — ).  Also  shown 
IS  the  curve  for  the  applied  osdilating  field  F^=  1. 5x1 0®V/m(-*-*-)  with  frequency  1.8x10” 
rad/s.  The  curve  (■— )  Is  a  normalized  version  of  the  curve  for  (-*-*-)• 


J  (10*A/rn) 


J  (10^A/rn) 
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Appendix  E 


100  by  100  pixels  can  take  several  minutes  to  display  as  a  three  dimensional  plot 
on  a  typicaJ  microcomputer.  The  same  plot  takes  30  seconds  on  this  system. 
If  twenty  images  are  recorded  in  a  single  day,  each  one  being  say,  500  by  500 
pixels,  displaying  that  data  on  a  simple  microcomputer  would  require  throwing 
away  every  other  pixel  and  it  would  still  take  hours  to  review  a  days  worth  of 
data. 

The  last  function  of  the  computer  is  control  of  the  STM  scan.  This  requires 
a  special  output  device  to  send  the  correct  electronic  signals  to  the  electronics 
controlling  the  STM,  There  are  a  number  of  possible  ways  to  do  this.  There  is 
an  instrument  bus,  IEEE  488,  or  direct  digital  I/O,  or  an  analog  control  signal 
generated  with  a  digital  to  analog  converter.  Each  method  has  advantages  and 
disadvantages.  The  D/A  method  was  chosen  as  a  quick  and  dirty  approach. 
It’s  cheap,  and  can  interface  to  simple  analog  circuitry  controlling  the  STM.  It’s 
disadvantage  is  noise.  At  some  point  in  the  future,  a  better  but  more  expensive 
option  may  be  used  in  it’s  place. 

The  requirements  for  speed  and  accuracy  are  roughly  the  same  as  those  for 
the  A/D  converter  used  to  collect  data.  Now  let’s  outline  the  actual  system 
built  to  meet  these  functions. 


3  Meeting  the  Requirements 

3.1  The  Outline 

This  section  will  explain  how  the  above  requirements  were  met.  The  overall 
computer  system  for  the  STM  is  shown  in  figure  1,  This  paper  will  give  the 
details  of  the  hardware  shown  in  figure  1. 

3.2  The  Chassis 

The  PC  tower  was  chosen  as  the  basis  for  this  system.  The  reasons  were  expan¬ 
sion  capacity,  availability  of  boards,  and  cost.  While  the  PC  bus  is  one  of  the 
slower  ones  available,  CPU  speeds  and  video  speeds  have  increased  enough  in 
the  last  couple  of  years,  to  allow  the  PC  to  be  the  basis  of  our  system.  Only  a 
couple  of  years  ago,  STM  would  have  required  a  much  more  expensive  bus  such 
as  VME.  The  tower  case  was  needed  to  allow  enough  cards  to  be  installed  (6 
expansion  cards)  and  to  supply  enough  power. 

3.3  The  CPU 

The  central  processor  unit  (CPU)  card  is  critical.  In  order  to  meet  the  speed 
requirement  outlined  above,  the  CPU  had  to  be  fast  and  able  to  address  large 
amounts  of  memory.  High  speed  floating  point  was  also  required.  The  CPU 
card  chosen  has  an  Intel  486DX  CPU,  running  at  66  MHz.  This  is  an  minimum. 
A  100  MHz  Pentium  would  be  better. 
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3A  Extra  Memory 

Experience  had  already  shown  that  4  megabyte  of  memory  is  not  enough.  In 
order  to  be  able  to  handle  low  resolution  images  8  megabyte  is  just  enough.  So 
8  megabyte  of  extra  memory  has  been  installed.  The  total  RAM  in  the  system 
is  now  20  megabutes.  Xwindows  is  the  largest  user  of  RAM.  If  the  system  were 
used  without  Xwindows,  a  system  with  8  megabytes  of  RAM  would  do. 

3.5  Disk  and  Tape 

There  are  four  disk  drives  connected  to  this  system.  Two  hard  disk  with  a 
capacity  of  over  700  megabytes  and  a  two  floppy  drives,  one  5  inch  and  one  3 
inch  drive.  The  speed  of  disk  access  is  important  when  large  volumes  of  data 
must  be  reviewed  or  processed.  The  current  system  uses  an  IDE  interface.  The 
system  could  be  improved  with  SCSI  disk.  The  floppy  drive  is  used  to  make 
backup  copies  of  data  once  it  is  no  longer  needed  on-line.  A  tape  drive  which 
connects  through  the  floppy  controller  is  used  for  major  system  backups. 

The  capacity  of  the  hard  disk  will  be  a  problem  if  high  resolution  data  taking 
begins.  A  100  by  100  pixel  image  only  requires  20000  bytes  of  storage.  A  single 
image  of  500  by  500  pixels  requires  1/2  megabyte  to  store.  That  is  one  floppy 
disk  per  image.  While  the  current  system  can  just  get  by,  a  larger  hard  disk  and 
a  SCSI  tape  backup  would  be  a  good  idea.  These  can  be  added  to  the  current 
system. 

The  base  system,  chassis,  mother  board,  disks  and  tape  drive  are  all  rather 
generic  parts. 

3.6  Video  card 

The  graphics  requirements  are  currently  being  met  by  a  Diamond  “Viper”  SV¬ 
GA  graphics  caud,  and  a  14  inch  color  multiscan  monitor.  Both  of  these  devices 
have  a  resolution  of  512  pixels.  The  graphics  board  has  it’s  own  memory  and 
processing  firmware.  The  video  card  was  a  mistake.  Diamond  has  a  closed 
architecture  amd  hence  is  not  supported  by  second  source  software.  It  should 
be  replaced  with  a  better  card  such  as  the  ATI  “Wonder”  card.  The  monitor 
should  also  be  upgraded  to  at  least  a  17  inch  monitor  to  take  full  advantage 
the  Xwindows  environment.  If  the  mother  board  is  upgraded,  a  PCI  bus  with 
ATI  video  card  is  recommended  for  best  results  using  Xwindows.  The  graph¬ 
ics  firmware  supports  at  least  256  colors  or  256  grey  levels,  this  is  important 
for  image  display.  It  also  requires  an  analog  RGB  signal  to  work  at  such  high 
resolution.  It  takes  a  special  monitor  to  use  these  signals. 

3.7  Analog  and  Digital 

There  are  both  an  analog  to  digital  and  digital  to  analog  boards  in  this  system. 
Both  of  these  boards  were  supplied  by  Computer  Boards  Inc.  [?].  The  anaJog 
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to  digital  board  is  a  12  bit,  50  KHz  device.  The  D/A  board  is  100  KHz  and 
resolution  16  bits. 

The  A/D  boeo’d  is  used  to  record  tunneling  current.  The  D/A  board  controls 
the  STM  scans  by  suppling  analog  X,  Y  and  Z  signals  to  a  high  voltage  amplifier 
which  drives  the  piezoelectric  tube.  Also  the  D/A  board  supplies  the  bias  voltage 
to  the  tip. 

These  boards  are  just  fast  enough  to  support  a  256  pixel  image  scan.  To 
push  to  higher  resolution  imaging,  A/D  board  will  need  be  replaced  at  some 
point  in  the  future  with  a  16  bit,  100  kHz  board.  The  12  bit  resolution  of  the 
A/D  board  is  a  problem.  Twelve  bits  gives  a  resolution  of  1  part  in  4096.  But 
that  assumes  you  make  use  of  the  total  range  of  the  board.  In  practice,  this 
is  not  the  case.  If  the  board  has  a  10  volt  range,  and  the  signals  from  a  given 
experiment  only  span  1  volt,  then  the  effective  resolution  of  such  a  board  is  only 
1  part  in  400.  This  is  stiU  ok  for  pictures,  but  it  can  add  unwanted  noise  to  the 
signaJ  if  statistical  analysis  is  to  be  performed  on  the  data. 

3.8  The  Outside  World 

The  current  system  has  two  serial  lines  running  to  the  outside  world.  One  line 
is  connected  serial  mouse  and  is  used  as  input  for  Xwindows.  The  second  line  is 
connected  to  a  standard  ASCI  terminal  and  can  be  used  to  log  into  the  system 
to  do  background  task. 

There  is  also  a  standard  ethernet  card  connecting  the  system  to  hypespace. 
The  system  allows  remote  users  to  log  in  from  anywhere  on  the  internet.  It  also 
supports  remote  Xtenninals. 

3.9  The  Software 

Of  course  the  real  work  in  getting  all  this  hardware  to  function  together  as  a 
system  is  in  the  software.  OS9000/PC  was  the  operating  system  of  choice.  It 
is  small,  fast,  powerful  and  relatively  bug  free.  It  supports  C,  FORTRAN  and 
486  assembly  as  the  programming  languages.  It  is  siniilar  to  the  UNIX  system 
but  is  much  more  compact.  It  also  supports  real  time  operation,  which  is  a  real 
plus  in  a  laboratory  environment. 

After  the  operating  system,  the  next  level  of  software  is  the  system  of  hard¬ 
ware  drivers.  Each  of  these  programs  is  required  to  interface  the  operating 
system  and  applications  programs  to  specific  devices  in  a  standardized  way.  S- 
ince  this  computer  system  was  put  together  with  parts  from  several  different 
vendors,  some  of  the  device  drivers  (D/A  and  A/D)  had  to  be  developed  here. 
Software  to  support  the  serial  lines,  ethernet  and  Xwindows  were  all  supplied 
by  Mlcroware. 

The  next  level  of  the  software  is  applications  program  libraries.  There  are  li¬ 
braries  to  support  devices  such  as  the  graphics  processor  and  mouse,  Xwindows. 
Internet  support  libraries  are  supplied  by  Microware. 
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Finally,  there  is  the  applications  software.  These  programs  include  various 
utilities,  image  processing  and  display,  data  acquisition  and  STM  control. 

The  serial  links  to  the  outside  world  have  been  functioning  for  some  time 
and  work  well.  Ethernet  access  includes  telnet  and  FTP.  The  graphics  system 
is  functioning  and  emulates  a  black  and  with  Xterminal.  Both  the  Xserver  and 
the  video  card  must  be  upgraded  to  display  color  images  on  the  system  monitor. 
Currently,  images  are  displayed  on  a  remote  HP  Xwindow  workstation.  There  is 
still  a  lot  of  work  to  do  at  the  applications  level  to  develop  better  displays.  De> 
vice  drivers  for  the  A/D  and  D/A  boards  have  been  developed.  Data  acquisition 
and  STM  control  programs  have  been  written  and  improvements  continue. 

All  of  this  software  must  fit  into  some  sort  of  unified  whole.  The  goal  is 
to  acquire  and  analyze  STM  data.  Since  higher  resolution  imaging  is  the  main 
thrust  of  this  effort,  speed  and  efficiency  of  both  hardware  and  software  are 
very  critical.  This  requires  that  the  software  effort  be  well  organized  and  the 
programs  must  be  well  written. 
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Figure 


A  scan  of  a  planar  gold  surface  taken  immediately  after 
forming  a  17  X  44  Angstrom  pit. 
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ABSTRACT 

The  spatial  and  temporal  temperature  distribution  for  a  conical  sh;q)e  metal 

under  laser  iiradiance  is  obtained  analytically  using  the  Green  function  method. 
In  this  study  the  tip  is  modeled  as  an  infinite  cone  of  half  angle  Oowhich  can  vary 

between  0  and  it/2.  The  full  three  dimensional  heat  diffusion  equation  is  solved 
simultaneously  with  the  Fourier  equation  for  the  heat  flux  assuming  no  radiation  losses. 
The  general  solution  is  obtained  for  an  arbitrary  temporal  and  spatial  distribution  of  the 
irradiance.  With  uniform  irradiance  the  temperature  rise  in  the  conical  tip  varies  almost 
linearly  with  the  laser  intensity  and  heating  time  and  d^ends  strongly  on  cone  angle 
and  thermal  properties.  For  a  tungsten  tip,  the  temperature  increase  is  about  2-3  orders 

of  magnitude  above  ambient  for  a  typical  laser  intensity  of  about  IMW/cm^  in  a  time 
of  a  few  hundred  nanoseconds.  To  study  Joule  heating  and  the  Thomson  thermodectric 
effect  in  the  STM  junction  we  first  obtained  the  current  distribution  inside  the  emitter. 
Using  this  as  a  volume  heat  source,  we  used  the  Green  function  method  to  solve  the 
heat  conduction  equati<Hi.  The  resulting  temperature  gradient  was  then  used  to  obtain 

the  thermoelectric  potential.  For  a  tungsten  tip  at "  1000  °  K,  the  thermoelectric  effect 
generates  a  bias  vdtage  cm  the  order  of  10  mV. 
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1.  INTRODUCTION 


Recent  experiments  on  laser  frequency  mixing  [1-3]  and  tunneling  time 
measunnents[4]  have  prompted  interest  in  thermal  effects  due  to  laser  irradiation  of 
an  STM  junction.  When  the  laser  field  interacts  with  an  STM.  the  tip  and  base  act  as  an 
antenna  system  fcx'  detection,  harmonic  generation  and  mixing  [5].  The  resulting 
nonlinear  1-V  characteristics  of  the  STM  diode  junction  produce  both  a  rectified  dc 
and  ac  tunneling  voltage  bias.  In  addition,  the  thermoelectric  effects  may  also  contribute 
to  the  bias  potentials  in  the  junction.  To  investigate  these  latter  effects  it  is  necessary  to 
determine  the  laser  induced  temperature  dbtribution  in  the  junction.  The  resulting 
temperature  gradient  inside  the  emitter  as  well  as  the  rectified  dc  current  through 
the  STM  junction  can  produce  both  resistive  (Joule)  heating  [6]  and  a  Thomson 
thermoelectric  effect  due  to  the  absorption  or  evolution  of  heat  Thus  the  STM  junction 
can  be  heated  or  cooled  by  the  reversible  Thomsrm  heat  depending  on  die  directiem  of 
tunneling  current  flow.  These  heating  or  cooling  effects  modify  the  temperature 
distribution  inside  the  emitter  and  anode  with  the  result  that  a  thermoelectric  potential 
can  be  generated  by  a  temperature  gradient  between  the  tip  and  base.  These  heating 
effects  also  produce  thermal  expansion  [7]  so  the  tunneling  gap  width  will  be  modulated 
influencing  the  tunneling  current  To  find  the  temperature  distribution  inside  the  tip 
(and  anode  for  0=5^)  due  to  Joule  and  Thomson  heat  we  solved  the  Laplace  and 
continuity  equations  to  obtain  the  current  distribution.  To  avoid  mathematical 
singularities  and  replicate  the  current  distributiem  for  a  realistic  tip,  the  infinite  cone 
model  is  truncated  at  some  apprqpriate  distance  from  the  apex.  The  temperature 
distribution  is  then  obtained  by  solving  the  heat  conduction  equation  using  the  (Jreen 
function  method  with  surface  and  volume  generation  of  heat  Knowing  the  temperature 
gradient  the  laser  induced  thermoelectric  potential  in  the  STM  junction  is  calculated. 

In  section  11,  we  outline  the  calculation  for  the  temperature  distributimi  in  an  infinite 
and  truncated  conical  emitter  using  the  (jreen  function.  The  solution  for  surface  and 
volume  generation  of  heat  by  laser  irradiation  and  by  resistive  and/or  Thomsem  beat 
respectively  are  also  given  in  this  sectiem .  The  thermoelectric  vedtage  generated  by  heat 


and  current  sources  is  presented  in  section  m  and  numerical  result  and  discussion  is 
given  in  section  IV. 

n.  THE  CALCULATION  OF  THE  TEMPERATURE  DISTRIBUTION 

In  our  analysb  we  assume  that  the  emitter  tip  is  an  infinite  metallic  cone  of  angle  Bg 

which  can  varies  0  to  Wl  [see  The  emitter  is  immersed  in  a  uniform  azimuthally 
symmetric  laser  beam  incident  normally  on  the  cone  surface.  We  here  only  briefly 
summarize  the  details  of  the  calculation  of  the  temperature  distribution  with  surface  and 
volume  sources  of  heat  A  detailed  description  of  this  analysis  b  to  be  publbhed 
elsewhere. 

To  obtain  the  temperature  dbtribution  T(r,t).  we  solve  the  inhomogeneous  heat 
conduction  equation 


F<r.t) 

K 


(1) 


where  r  and  t  are  spatial  and  time  cocMrdinates,  k  b  the  thermal  diffusivity,  K  b  the 
thermal  conductivity  (k  and  K  are  assumed  as  constant)  and  F(r,t)  b  the  rate  of  heat 
supplied  per  unit  volume  per  unit  time. 

We  use  the  Green  fimction  method  to  obtain  a  general  solution  of  eq.  <1). 

The  procedure  b ; 

1.  Solve  for  the  Green  functicm  G<r,Ur',0  f<v  a  point  source  with  N^unann  boundary 
conditions 


n.  VG<r.Lr',O=0  (2) 

where  prime  denotes  the  source,  and  n  b  an  outward  pointing  normal  to  the  surface. 
The  Green  function  satisfies 
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V^(Kr,ur't>- - - — - 


.  4k  6(r^  6(1-0 


(3) 


It  is  assumed  that  both  volume  and  surface  generation  are  possible. 

2.  Apply  Green's  theorem  and  the  principle  of  superposition  to  obtain  the  solution  for 
an  arbitrary  source  distribution. 

For  volume  and  surface  generation  the  resulting  temperature  distribution  is  given  by  [8] 


T(r.t)=P^dtJ‘^G<r.tj'.0^^^  k  dS  +JJ^dt'JyG<r,ty,t'> 


‘''^^kdV  (4) 


The  boundary  condition  to  be  satisfied  fcv  surface  generation  is 
aT(r'.t')  Q<r',t') 

an'  K  ^  ^ 

where,  in  our  case,  CKr'tO  is  the  heat  produced  per  unit  time  per  unit  area  due  to  laser 
irradiation  at  the  surface. 

It  is  assumed  that  in  general  the  laser  intensity  has  Gaussian  spatial  and  temporal 
distribution: 

Q(r'.0=Qo<l-RO^))  exp<-<r'/  D^)  eKp<-<t'/ri>)S  (6) 

where  Qo  b  the  maximum  intensity  of  the  laser  beam.  R(X)  is  the  reflectivity  and  Tq,  D 

are  temporal  and  spatial  pulse  widths.  The  reflectivity  is  a  function  of  the  wavelength  X. 

We  first  consider  the  emitter  to  be  modelded  as  an  infinite  cone.  Subsequently 
we  will  concern  ourselves  only  with  the  temperature  distribution  near  the  dpex  of 
the  emitter.  It  can  be  shown  that  the  (rreen  function  for  this  case  is  [9]: 

G<r.ty  ,0=  -  ,  I  cxp[  V+r' Vk(t-t')] 

4ic  kvt-i  v«  /  n 


^  *0+1/2  <rr'/2k(t-t')) 


dnd^ 


<7) 


»»-l*o 


*0+1/2  ntiodifled  Bessel  fu0ctio0  aod  PdQI)  u  the  Legeodre  fu0ctio0. 

.dP„0i)| 


Tbe  summatio0  over  0  is  for  succesive  roots  of - 


d|i 


0  greater  tha0  -1/2. 


i»*m> 


For  the  fmite  (i.e..  truocated)  co0e  case  [fig.2],  we  assume,  as  m  tbe  mfioite  co0e 
case,  a  uoifcxm  azimuthally  synuoetric  laser  irradiatio0  is  i0cide0t  00  the  ccme  surface. 
The  cone  leugth  is  assumed  large  euough  compared  to  the  real  emitter  radius  so  that 
the  upper  part  of  cooe  (at  radius  r=d  io  fig.2)  is  kept  at  a  coostaot  teniperature  hy  a  heat 
reservoir.  The  result  b 


(Kr,ty,t>-EE 

um 


C^'^1)  Jn*1  ^Ynm  f)lD4-l/2^Ynni  ^ 

d^tWY.md)f 


-  cxp(.tL.«.-0) 


90  dp 


(8) 


where  J  ,  „(x)  b  a  Bessel  limctioo  and  tbe  siunmatio0  over  m  b  for  roots  of 
n+1/2' 

■*n+l/2  ^Tnm  ®' 

n.l.  SURFACE  GENERATION  OF  HEAT 


For  a  U0iform  q;>atial  and  temporal  dhtributicm  of  CW  laser  radbtion  at  the  surface 
of  an  infinite  cone  model  the  temperature  dbtribution  [9]  b 


T(r,e,t)=EE  Qoa-R<X)) 

na 


qn+l)  rIXa+0/2+l) 

4K  IXa+1)  IXa+n+3/2) 
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(9) 


?  PnOi) 

r(a+n/2-l/2  r— ) - 

4ln 


where  IXx) » the  complete  Gamma  function  and  r(x;y)  >*  the  incomplete  gamma 
function. 

In  the  case  for  which  the  cone  degrades  into  a  plane  (0o«x/2)  the  temperature 
distribution  along  the  axis  (6=0)  is 


T<r,t)= 


Ood-RO.)) 


rIX-l/2;^> 


4K 

Qo<l-R<X)) 

4K 


1X3/2) 


(10) 

(11) 


where  ierfc  (x)=J*«fc  (z)  dz  and  erfc  is  the  errcx  function.  The  planar  case 
(xxrespond  to  the  temperature  distribiution  in  the  laser  inadiated  anode. 

n2  THE  VOLUME  GENERATION  OF  HEAT 

n2.1  THE  CALCULATION  OF  CURRENT  DISTRIBUTION 

In  this  section  we  will  consider  the  temperature  distributkai  in  the  emitter  tip  by 
resistive  generation  of  heat  due  to  a  tunneling  current  flow  across  the  STM  junction. 

To  avoid  a  singularity  in  the  cunent  density  at  the  zpex  of  the  cone,  we  assume  that  in 

the  region  denoted  asl  for0sr2sc(in  fig.2),  the  constant  current  density  - 10^- 10* 

amp/cm^  flows  into  the  r^cm  through  the  finite  conical  t^  surfac^  This  ;q}proximation 
is  not  unreasonable  because  it  assume  a  more  realistic  current  distribution  over  a  finite 
area  rather  through  a  singular  point  Another  justification  for  fiiis  tqiproximation  is 

that  the  field  in  this  r^ion  decays  so  rs^idly  that  for  r  >  c.  it  becomes  negligible  and 
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Pni<l*')Po2<t*') 


902  9n 


IM-Mo 


(13) 


n^.3  THOMSON  HEAT  AS  A  SOURCE  TERM 


When  both  an  the  electric  and  heat  current  exist  at  the  same  time  in  the  volume  of 
a  conductor,  there  can  also  be  Thomson  heat  produced  or  absorbed  in  addition  to 
resistative  or  Joule  beat  The  Thomson  heat  evolved  (or  absorbed)  per  unit  volume 
per  per  unit  time  is  given  J«VT  (When  the  Thomson  beat  is  absorbed,  the  Tbomscm 

coefficient  is  positive  if  J  and  VT  are  parallel)  [11].  As  a  source  of  volume  generation 

of  heat  we  use  the  current  density  obtained  in  sec.  II.2.1  and  the  temperature  distributimi 
given  in  sec.  n.l.  The  source  distribution  (le.,  the  Thomson  heat  evolved  in  this  system) 
is  F(r,t)=o^  J«VT.  The  temperature  distribution,  with  Thomson  beat  as  a  source,  can  be 

found  from  eq.  (4): 


Tr,«,w,  I  IXa+ng^l) 

^  4K  Q  r(<i+l)  r(a+D+3/2) 


JyJ^  G(r,ty,0 


X  r(a+n/2-l/2; - ) 

4k(t-0 


<1-1*^ 


fcdV.n  Qod-Ra))!  _  -  _ 

z  K  d  r  +  o_ - I  £ 

3^n(l*)l  ^  8*  K(l-|iJ)^  n^iija 


dndp 


x(2ni+lX2n2+l> 


IXct+nj/Z+l) 


IXo+l)IXa+ni+3/2> 


Pni<l‘')Pn2<l*') 

X  IXa+11,/2-1/2; - ) - 

^  4k(t-0 


J  dtj*^  CKr.tj',0^ 


kdV 


(n2  c“2*^+d^\n2+l)c'^ 
(14) 


dni  9|i 


l»-m  ,  dn2d|i 


-8- 


no  current  is  assumed  to  flow  through  the  surface  of  the  cone  in  the  region  2  (see  flg.2) 

[10].  The  spatial  polar  coordinate  3xl0'^cm  defines  the  lower  boundary  of  region  2. 

We  also  truncate  the  cone  at  r^d"  1mm.  To  obtain  the  current  density  distribution  in 
region  2,  we  solve  the  Laplace  equation  for  the  potential  of  the  current  field  J=>o  W 
satisfying  the  continuity  equation  V*J=0.  The  resulting  current  density  distribution  in 
region  2  (c  s  r  s  d)  is  given  by 


J2(r.0)=-oW 

(2n-H)  (nr°'^+d^'^\n-H)r'°'^  Pn(H) 


=z 


(12) 


dn  dp 


where  V  is  the  electric  potential,  o  is  the  electrical  conductivity  and  the  current 
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IsJqX  c  sinOo  in  this  regkm  is  constant 


n.2.2  JOULE  HEAT  AS  A  SOURCE  TERM 


In  this  section  we  calculate  the  temperature  distribution  assuming  the  Joule  heat 

as  a  source  tenn.  This  vc^ume  source  term  is  F(r,t>=  J^p,  where  p  is  the  electrical 
resistivity.  The  temperature  distribution  inside  the  cone  emitter  is  obtain  from  eq.(4): 


T(r,0.t)=J^dt'J^^G(r.tr'.O^kdV+£  I  (2n^+lX2n2+l> 


I^P 


fdt'f  G(r.tr'.t'> 
^  0 


(nic**l*^4d^^*^(nj^+l)c  °1)  (n2  c°2*^+d^^^(n2+l)c*°2) 


In  the  next  section,  we  describe  the  calculation  of  the  thermoelectric  voltage  due  to 
both  surface  and  volume  sources  of  heat 

ra.  THERMOELECTRIC  POTENTIAL 

We  have  previously  considered  the  ^atial  and  temporal  temperature  distribution 
due  to  the  laser  irradiance  of  the  STM  junction.  This  produces  not  only  a  temperature 
gradiant  due  to  surface  generation  of  heat  but  also  a  rectified  dc  cunent  component 
through  the  junction.  The  thermal  energy  in  the  tip  consists  of  beat  brought  into 
the  volume  by  the  beat  conduction  through  the  surface  as  well  as  resistive  heat  due 
to  the  current  density  established  by  the  laser  field  even  in  absence  of  an  external  bias 
voltage.  In  addition,  there  b  a  Thomson  beat  due  to  the  temperature  gradient  as  well 
as  current  flow  in  the  junction.  Thb  heat  can  be  either  absorbed  or  evolved  into  the 
system  influeix;ing  the  temperature  gradient  across  the  juctiotL  The  thermoelectric 
potential  generated  by  the  field  due  to  the  temperature  gradiant  between  two  poinb  in 
the  presence  of  current  can  be  written  as  [11] 


(15) 

-/OjV,T* 

(16) 

.  II 

<17) 

=TiSai)-ToS<To)-/^^  SdT 
^0 

(18) 

dS 

where  o_  b  related  to  the  thermoelectric  power  S  by  o_=  T— . 

*  *  dT 

Therefore  the  thermoelectric  potential  produced  across  the  junction  b 
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Vj;,K  -TiSeCTi)  -  T2S,<T2)-  FlS.  dT  S.  dT 

*0  Tq 


(19) 


where  T^and  T2  are  the  temperatures  of  the  emitter  and  anode  base,  and 

are  the  absolute  thermoelectric  powers  of  the  emitter  and  anode  base,  respectively. 

It  is  assumed  that  at  suMciently  large  distances  from  the  apex. 

For  the  sake  of  completness  we  also  consider  the  current  density  due  to  the  ^atial 
dependence  of  the  chemical  potential  Since  this  quantity  is  a  function  of  the 
temperature  the  difrierential  variation  in  chemical  potential  with  temperature  in  the 
junction  can  also  produce  a  thermoelectric  voltage  in  the  juction  [12].  Therefore  the  total 
field  produced  in  the  system  is 


GO) 


wb^e  J/o  is  the  field  produced  by  the  laser  irradiation  in  the  absence  of  an  external 
applied  voltage.  The  thermoelectric  potential  generated  due  to  the  effective 

field  Ecg=-Vf^-OxVfT  in  the  STM  junction  is 
c 


GD 


or 


V  —  dT  -  ^  dT-f ^2  T  ^  dT) 

c^TodT  ^TodT  dT  dT  ^ 


(22) 

wbae  ^  and  ^  are  the  chemical  potentials  of  the  emitter  and  anode,  respectively. 
In  the  case  where  T]^=:T2*  T  (T  is  the  local  equilibrium  temperature  at  the  junction 
and  assuming  To=0)  the  effective  thermoelectric  potential  becomes 
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“  ■  e  ^  T[Sla>S2C^)^  /^(Se  -S,  )  dT  <Z3) 

However,  for  the  present  application,  this  contribution  is  of  second  •order  and  will 
not  be  included  in  the  determination  of 

IV.  DISCUSSION  AND  CONCLUSIONS 

In  this  pq;>er  we  have  studied  the  thermal  effects  in  an  STM  junction  under 
laser  inadiance.  For  a  tungsten  emitter  modeled  as  a  cone,  we  have  used  a  Green 
fimction  method  to  obtain  analytic  expressions  for  the  temperature  distribution  with 
both  surface  and  volume  sources  of  beat  Using  the  laser  irradiance,  resistive  and 
Thomson  beat  as  the  surface  and  volume  sources  respectively,  we  have  calculated  the 
thermoelectric  potential  amtribution  to  the  bias  voltage  across  the  junction.  Some 
q;>ecific  results  we  have  obtained  are  sununarized  below. 

The  tenq>erature  rise  in  the  conical  ti^  doe  to  surface  generation  of  beat  varies 
almost  linearly  with  laser  intensity  and  heating  time  and  depends  strongly  on  the  cone 
an^e  and  thermal  parameters  [fig.3] .  In  a  beating  time  of  about  120  ns  with  a  CW  laser 

intensity  of  1  MW/cm^ ,  the  temperature  rise  of  a  conical  emitter  with  half-  cone  angle 

of  15^  is  about  1100  ^  K  near  the  zpex  (within  1  nm)  and  drops  r^idly  for  distance 
beyond  a  micron  from  the  ^x  [fig.4]. 

The  temperature  distribution  due  to  Joule  heat  produced  by  the  laser-  induced 
tunneling  current  dqiends  strongly  on  the  magnitude  of  die  induced  current  density 

[fig.5].  For  densities  in  the  range  lOl^-lO^  aaap/ca?  the  temperature  increase  is 
negjegible  because  the  current  density  is  so  small  and  the  thermal  conductivity  of 

the  material  is  high.  When  the  current  density  is  about  10^  amp/cm^  ,  the  temperature 

increase  is  about  50  ®  K.  This  temperature  change  is  iq>ivoximately  uniform  near  the  tqi 
apex  and  then  begin  to  drop  riqiidly  for  distance  beyond  100  nm.  When  die  current 

density  is  greater  than  5  x  10^  amp/cm^  the  temperature  of  the  conical  emitter  (half  cone 
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angle  of  IS  increasea  dramatically  and  is  on  the  order  of  2000  These  temperatures 
already  begin  to  exceed  the  melting  temperature  of  most  metals.  For  comparision,  Tn,  for 

tungsten  is  3683  °  K. 

The  temperature  change  due  to  Thomson  heat  generated  or  absorbed  by  the 
system  is  relatively  small  compared  to  the  both  the  laser  induced  surface  heating  and 
the  resistive  heating  due  to  the  tunneling  current  because  the  current  flow  is  small  and 
consequently  the  heat  transport  through  the  emitter  is  so  small.  Therfore  the  heating  or 
cooling  effect  due  to  the  Thomson  heats  is  negligible  compare  to  the  surface  and  J(Mile 
heating  effects. 

Lastly  the  thermoelectric  potential  is  due  princ4»ally  to  surface  heating  and  depends 
on  the  intensity  of  the  laser  irradiation.  The  induced  thermoelectric  voltage  due  to 

surface  generation  of  heat  for  a  laser  intensity  of  12  MW/cm^  is  about  10  mV  [fig.6]. 
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Figure  Cq)tioos 


Fig.l  The  infinite  cone  model.  The  half-  cone  angle  6o  varies  from  0®  to  90®. 

Fig.2  The  flinite  cone  model.  The  cone  is  truncated  at  r=d  (- 1  mm).  The  active  region 

for  electron  emission  is  the  conical  surface  between  r=0  and  r=c  ("Sx  10  ^  cm). 

Fig.3  The  laser-  induced  tempereture  rise  near  the  apex  (10  A®)  vs.  heating  time  for 
the  infinite  cone  model.  The  solid  line  b  for  a  laser  intensity  of  1.6MW/cm  and  the 
dashed  line  b  for  1  MW/cm^.  The  cone  half  angle  b:(a)15®  (b)30®  (c)90®.  The  values  of 
k=0.68  cm2  8*^,  k=1.78  W/cm®K,  and  R=0.68. 

Fig.4  The  laser-induced  temperature  rise  vs.  radial  dbtance  along  the  emitter  axb 
for  different  cone  angles  after  a  heating  time  of  120  ns  with  a  laser  irradiance  of 

1  MW/cm2.  The  same  value  of  thermal  parameters  as  in  fig.2  are  used. 

Fig.S  The  temperature  rise  due  to  Joule  heating  vs.  radial  dbtance  along  the  axb 
of  the  emitter  for  different  current  densities.  The  heating  time  b  100  ns.  The  resbtivity 

of  the  tungsten  sample  was  taken  to  be  5x  10  ^  ohm  cm. 

Fig.6  The  laser  induced  thermoelectric  voltage  vs.  energy  flux  of  CW  laser  beam  with 
different  cone  angle.  The  surface  generation  of  heating  effect  b  only  considered  in  thb 
data  and  the  heating  time  b  100  ns. 
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